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ABSTRACT 

Multi-star systems are common, yet little is known about a stellar companion’s influence on the 
formation and evolution of planetary systems. For instance, stellar companions may have facilitated 
the inward migration of hot Jupiters towards to their present day positions. Many observed short 
period gas giant planets also have orbits that are misaligned with respect to their star’s spin axis, 
which has also been attributed to the presence of a massive outer companion on a non-coplanar orbit. 

We present the results of a multi-band direct imaging survey using Keck NIRC2 to measure the 
fraction of short period gas giant planets found in multi-star systems. Over three years, we completed 
a survey of 50 targets (“Friends of Hot Jupiters”) with 27 targets showing some signature of multi¬ 
body interaction (misaligned or eccentric orbits) and 23 targets in a control sample (well-aligned and 
circular orbits). We report the masses, projected separations, and confirmed common proper motion 
for the 19 stellar companions found around 17 stars. Correcting for survey incompleteness, we report 
companion fractions of 48% ± 9%, 47% ± 12%, and 51% ± 13% in our total, misaligned/eccentric, and 
control samples, respectively. This total stellar companion fraction is 2.8 a larger than the fraction 
of field stars with companions approximately 50 — 2000 AU. We observe no correlation between 
misaligned/eccentric hot Jupiter systems and the incidence of stellar companions. Combining this 
result with our previous radial velocity survey, we determine that 72% ± 16% of hot Jupiters are part 
of multi-planet and/or multi-star systems. 

Subject headings: binaries: close — binaries: eclipsing - methods: observational planetary systems 
— planets and satellites: dynamical evolution and stability — techniques: high 
angular resolution 


1. INTRODUCTION 

Surveys of Solar-like stars within 25 pc indicate that 
appro ximately 44% are fou nd in multiple star sys¬ 
tems dRaghavan et al.1 l2010f) . At the same time, re¬ 
cent surveys have sought to quantify pla net occurrence 
rates a round solar-type F GK stars (e.g. iHoward et al.l 
l2012bl iFressin et aUljQljl ). However, the effects of ad¬ 
ditional stellar companions on the formation and subse¬ 
quent evolution of planetary systems are not well un¬ 
derstood. A stellar companion might di srupt p l anet 
form ation by stirring up the di sk (e.g. iMaver et al.1 


2005K . truncating the disk (e.g. iPichardo et al.l 12005 


Kraus eUdjUgOlg) or ejecting planets (e.g iKaib et, al I 


20131 iZuckermanl I2014T ) . Numerical simulations often 


fail to produce planets in binary star systems (e.g. 
Pichardo et al.1 120051 IMaver et al.l [20051 iThebault et al.l 


2006HFragner et ;dll2()l 111 suggesting that a stellar com- 
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panion can indeed hinder planet formation. On the other 
hand, analytic calculations predict that stellar compan- 
ions would have litt l e effect on plane tesimal growth (e.g. 


Batvgin et al.l 1201 It iRafikov 2013a,bT) and current sur- 
veys have found a number of planets in binary star 
systems (e.g. lEggenberger et al.l 120071: iRaghavan et al. 

201(1 IKaib et all 

20111: Orosz et al.l 2012a|]bi). In addi- 

tion, a stellar co 
grate via three-b 
Lidov mechanism 
ing in very small 

mpanion might cause planets to mi- 
ody interactions, such as the Kozai- 
or via other secular interactions, result- 

orbital distances (e.g. Malmberg et al. 

2007a| 

btlFabrvckv & Tremaine! 20071: IMorton & Johnson 

2014 

Naoz et al. 

2012L 20131: Tevssandier et al. 20131: 


Petrovichll2015t iStorch et al.ll2014l b However, the Kozai- 


Lidov mechanism is suppressed in multi-planet sys¬ 
tems because planet-planet interactions tend to prohibit 
the libration of the argument of peri helion characteris¬ 
tic of the Kozai-Li dov resonance (|Wu fe Murravl [20031 : 
iBatvgin et al.ll20lTI) . Finally, stellar companions can also 
bias our estimates of the properties of transiting planet 
systems by diluting the measured transit depth, resulting 
in an underestimate of the planet’s radius and a corre¬ 
sponding overestimate of its density. 

In this study we focus on a class of short-period gas 
giant planets known as “hot Jupiters.” These plan¬ 
ets could not have formed at their current locations, 
but must have migrated in from beyond the ice lines 
of their natal disks (e.g. iLin et all 1199 6). However, 
the mechanism(s) responsible for hot Jupiter migra¬ 
tion remain controversial. Current migration models in- 
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elude disk interactions (e.g. iGoldreich fc Tremainelll98d : 
iTanaka et al.l l2002t iLin fc Papaloizoul 1198611 and grav¬ 
itational interac tion with a third body, such as an¬ 
other planet (e.g.lChatt eriee et, al. 2008; Naeasaw a et, al l 


ottier planet le.g.ltm attcnc c et al.ll/UU S: iNasasaw a et, all 
gOOglJw^fcJnt hwic d 120111 iBeauae fc Nesvornvl l2012t 


Lithwick fc Wnl 12014 1 or a stellar companion. In gen¬ 


eral, isolated simple disk migration models produce hot 
Jupiters on circular orbits that are well aligned with the 
primary star’s spin axis, while migration due to a third 
body leads to hot Jupiters that are often eccentric and/or 
misaligned with the primary’s star spin axis. __ 

Surveys from the p ast few years (e.g. lWinn et al.l[2010l 
lAlbrecht et a,l.|[2012f) indicate misaligned hot Jupiters are 
common-18 out of the 53 hot Jupiters surveyed to date 
have obliquities that are inconsistent with zero at the 
three sigma level or higher. As a result, it has been 
argued that a significant fraction of hot Jupiters may 
have migrated via three-body interactions such as the the 
Kozai-Lidov effect, whi ch natura lly resu lts in large or- 
bital inclinations (e.g. iMorton fc Johnsonl 120111 : iLi et al.l 
(Ml). If stellar tides can bring misaligned hot Jupi ters 
back into alignment with the star’s spin axis (jDawsonl 
Ml), this fraction may be e ven higher than the cur- 
rent rate suggests. Conversely, iDawson et al.l ( 2015 b ar¬ 
gue that the lack of high eccentricity Jupiters at inter¬ 
mediate periods in the overall Kepler sample places a 
strict upper limit on the fraction of hot Jupiters that 
might have migrated via three-body processes. Mis¬ 
aligned hot Jupiters may also result from migration in 
a tilted disk, which could be caused by to rque from 
a distant stellar companion dBatvginl l2012f) . More¬ 
over, significant star-disk misalignments may naturally 
arise from the physical evolu tion of the star and the 
disk in a pertu rbe d syste m ((Batygin fc Adamsl 120131 ; 
iSnalding fc Batvginl l2014ah . This suggests that a hot 
Jupiter’s obliquity, which c an be measured with the 
Rossiter-McLauglilin effec t (IWinn et afl 120051) or via 
Doppler tomograp hy (e.g. iCollier Cameron et al.1 120101 : 
iBrown et al.l 12012) . might provide a clue to whether or 
not a third body has influenced the planetary system. 

Alternatively, planet-planet scattering could produce 
misaligned hot Jupit ers wi tho ut req uirin g the presen ce of 
a stellar companion. iDawson fc Murrav-Clavl (12013d also 
find evidence that high-eccentricity proto-hot Jupiters 
are more common around metal-rich stars, which pre¬ 
sumably are more likely to have formed multiple gas 
giant planets. Other studies have also suggested that 
protoplanetary disks in isolation might in fa ct be tilted 
by the chaotic nature of star formation (Bate ct al 


201C ), the primary star’s magnetic torques (lLai et al 


201 Id . and stellar surface modul ation by internal gravity 


waves ( Rogers et al.ll2012l 12013d . 

If a significant fraction of hot Jupiters migrate inward 
and acquire spin-orbit misalignments via three-body in¬ 
teractions, then this necessarily requires the presence of 
a massive outer planetary or stellar companion in these 
systems. However, there have not been any studies pub¬ 
lished to date that have provided a well-constrained es¬ 
timate of the frequency of bound stellar companions in 
hot Jupiter systems. A few stellar companions to transit¬ 
ing planet host stars were discovered serendipitously as 
part of studies intended to bette r characterize the tran¬ 
siting planet and its ho s t star ( e.g.lColl ier Cam eron et all 
120071 (Crossfield et al.1 l2012t iSing et al.1 12013T) ! Some 


other works, such as iNarita et al.l (l2010bl l2012d report 
directly imaged stellar companions from adaptive optics 
(AO) follow-up of known planets but only for one or two 
transiting gas giant planetary systems. The first system¬ 
atic surveys for stellar companions to transiting planet 
systems used the “Lucky imaging” method. These stud¬ 
ies focused exclusively on transiting hot Jupiter systems 
and th eir sa mple sizes w e re sm all: 14 in [Daemge^^^^' 


(2000), 16 in lFaedi et al.1 (|2013l) and 21 n lBergfors et a l 


( 2013 ). In addition, many of these surveys observed over¬ 
lapping target lists. 

More recently, there have been a series of studies fo¬ 
cusing on the sample of Kepler transit ing p l anet can¬ 
didate host stars. The two surveys by iLillo-Box et al.1 
(Ml Ml . covering 174 Kepler planet host candidates, 
is the largest “Lucky Imaging” search to date. Cur¬ 
rent state-of-the-art direct imaging surveys use adap¬ 
tive optics to achieve diffraction-limited im aging to allo w 
for be t ter detection and su r vey efficiency. lAdams et al.l 
(12012D. lAdams et al.l (12013D . iDressing et alT ( 20141) . and 
I Wang et al.l ( 2014bD obtained infrared AO images of 90, 
12, 87 , and 56 Kepler pl anet candidate hosts, respec¬ 
tively. lAdams et al. (|20 13l) also searc hed around 15 tran¬ 
siting planet hosts. lLaw et al.l (|2Qll recently published 
the first part of an optical campaign to search for com¬ 
panions around all Kepler planet candidate hosts using 
the Robo- AO instrument, with a n initial sample size of 
715 stars. lOilliland et al. (120151 ) searched for compan¬ 
ions around 23 Kepler planet candidate hosts using op¬ 
tical images from the Hubbl e Sp a ce Te lescope WFC3 in¬ 
strument. Finally, iHorch et al.l (I2014D used differential 
speckle imaging in two optical bandpasses to search for 
companions around 623 Kepler planet candidate hosts. 

Unlike previous surveys of hot Jupiters detected by 
ground-based transit surveys, these imaging surveys were 
intended to confirm the planetary nature of the tran¬ 
sits detected by Kepler and to correct the transit light 
curves for any dilution due to nearby stars, therefore 
ensuring accurate planetary radius estimates. Because 
the typical proper motions of the Kepler stars are quite 
small, these studies report relative brightness and pro¬ 
jected separation for companions but do not attempt to 
determine whether or not they are bound companions 
or background objects. The planetary systems in these 
surveys have a size distribution that reflects that of the 
Kepler survey as a whole, with the majority of systems 
consisting of sub-Neptune-sized transiting planets. 

In this work, we present a diffraction-limited direct 
imaging survey of close-in transiting gas giant planets or¬ 
biting bright, nearby stars, as part of the “Friends of Hot 
Jupiters” campaign. These systems are among the most 
favorable targets for the Rossiter-McLaughlin technique, 
and the majority of our targets have published measure¬ 
ments of their spin-orbit alignment. By focusing on this 
sample, we can directly test current hot Jupiter migra¬ 
tion models and investigate the origin of their observed 
spin-orbit misalignments by searching for massive, dis¬ 
tant companions in these systems. Our survey uses mul¬ 
tiple bandpasses and repeated observations spanning a 
several year baseline in order to determine whether any 
directly imaged companions are physically bound. We 
also use these same data to estimate companion masses 
and projected physical separations, which are required 
in order to evaluate the likelihood of specific dynamical 
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evolution scenarios for these systems. 

The Friends of Hot Jupiters survey uses multiple com¬ 
panion detection modes to search for planetary and stel¬ 
lar companions around exoplanetary systems. Our sam¬ 
ple consists of fifty-one exoplanetary systems that are 
known to host a transiting gas giant planet with masses 
between 0.06—11 Mj up and periods between 0.7—11 days. 
We divide this sample into two sub-samples, consisting 
of planets that are on misaligned and/or eccentric orbits 
and a control sample of planets on appare ntly circular, 
well-aligned orbits (see lKnutson et alJl20l4 for a full de¬ 
scription of the sample selection for this survey). We 
consider targets to be “misaligned” if they host planets 
with an eccentricity or spin-orbit alignm ent more than 
3 stan dard deviations away from zero. In lKnutson et ahl 
(120141) , we presented our search for long-term radial ve¬ 
locity (RV) accelerations due to distant massive plane¬ 
tary or stellar companions in these systems. We found a 
total companion occurrence rate of 51% ± 10% for com¬ 
panions with masses between 1 — 13 Mj up and orbital 
semimajor axes between 1 — 20 AU, with no evidence for 
a higher frequency of radial velocity companions in sys¬ 
tems with eccentric and/or misaligned gas giant planets. 
In a future paper we will present the results of a com¬ 
plementary search for close-in stellar companions using 
high resolution K band spectroscopy, which is primarily 
sensitive to K and M stars within OA^of the primary. 

In this paper, we present the results of our diffraction- 
limited direct imaging search. In Section [2j we describe 
our observations. In Section [3J we summarize the point 
spread function (PSF) fitting method used to calculate 
the brightness ratio and positions of the candidate stel¬ 
lar companions, as well as upper limits for compan¬ 
ions in systems with non-detections. We then determine 
whether or not the candidate companions share common 
proper motion with the primary, and estimate their pro¬ 
jected physical separations and masses. In Section [2 we 
discuss each system individually. In Section 0 we com¬ 
pare our estimated frequency for stellar companions to 
the results from previous surveys of planet-hosting and 
field star samples. Finally, in Section 1 we summarize 
our findings and discuss the implications of our measured 
companion fraction for the formation of hot Jupiter sys¬ 
tems. 


2. OBSERVATIONS 

During the AO phase of our survey, we collected data 
for 50 out of 51 FHJ systems with the NIRC2 instru¬ 
ment (instrument PI: Keith Matthews) on Keck II us¬ 
ing K band natural guide star adaptive optics imaging. 
We were not able to image one target, WASP-19, be¬ 
cause its declination, —45.7°, was too far south to ob¬ 
serve with Keck AO. Two of our target systems, HAT-P- 
8 and WASP-12, turned o ut to be triple system s, which 
we previously reported in lBechter et al.l (120141) . We ob¬ 
tained observations between February 2012 and October 
2014 and our observations are summarized in Table [1] 
We used the full array (1024x1024 pixel field of view) on 
the narrow camera setting (10 mas pixel” 1 ) to maximize 
our spatial resolution. However, for several bright tar¬ 
gets (as noted in Table [1]), we used a subarray to reduce 
integration times and avoid saturation. We used a three- 
point dither pattern to reduce the effects of the NIRC2 
array’s noisier lower left quadrant and instrumental noise 


levels while also preserving our sensitivity to companions 
with higher spatial separations. We aimed for a total of 
two minutes of on-target integration time per system in 
position angle mode, where the orientation of the image 
is kept constant on the detector as the telescope tracks. 
This technique allows us to detect companions with A K 
of approximately 8 at separations of approximately 1". 
For targets where a potential companion object is seen, 
we repeat the observations in at least one other filter, 
such as J or H, in order to obtain color information. We 
also follow up on targets with detected companions ap¬ 
proximately one or more years later to obtain K band 
astrometric measurements necessary to confirm that the 
companion is gravitationally bound via a common proper 
motion analysis. We elect to use K band rather than J 
or H band for our astrometry because the AO correction 
is superior at longer wavelengths. 

We calibrate our images using dome flat fields and dark 
frames. We also find and remove image artifacts. We flag 
flat field pixels that are less than 0.1 times the median 
as dead pixels and dark frame pixels that are more than 
lOer from the median as hot pixels. For each frame, we 
identify the remaining bad pixels as those with counts 
that are 8cr outliers compared to the counts in pixels in 
the surrounding 5x5 box. We replace all the flagged pix¬ 
els’ value with the median of the 5x5 box centered on the 
flagged pixel. We use these calibrated individual frames 
in all of our photometric and astrometric analyses. We 
limit our integrations to stay just below the nonlinear 
regime for the NIRC2 detector, and use Poisson statis¬ 
tics to determine the uncertainty in our counts. We also 
create a single, reduced image by aligning the individual 
frames so that the target star is in the same position and 
then combine using a median stack. We use the stacked 
image for our sensitivity calculations. 


3. ANALYSIS OF COMPANION PROPERTIES 
3.1. Detections 

We find 15 binary systems and 2 triple systems, for a 
total of 17 multi-star systems, out of the 50 systems with 
AO observations. We show one median-stacked K-band 
image for each of these detections in Figure [Q Table [2] 
summarizes the stellar parameters for all Friends of Hot 
Jupiters survey targets and the number of companions 
found around each star. 

To measure the flux ratio and on-sky separations 
for each system, we fit a multiple-source point spread 
fun ction (PSF) t o eac h calibrated frame. Follow¬ 
ing [Bechter^eFjil] (12014D . we choose to model the PSF 
with as a Moffat function with a Gaussian component, 






+ 7 i exp 



( 1 ) 


where N* is the number of stars in the image (either 2 
or 3); r, = yjxf + yf is the distance from the i -th star; 

are parameters that vary with each star and 
determine the position of the star and the amplitude of 
the PSF; /3 takes a single value for all stars and sets 
the exponent of the Moffat contribution; r s and w each 
take a single value for all stars and determine the width 
of the Moffat and Gaussian portions of the PSF, respec¬ 
tively; and b is the background sky level. Hence, the total 
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Figure 1 . Median-stacked K-band images showing the Friends of Hot Jupiters survey targets with detected and confirmed stellar com¬ 
panions (marked with arrows). The first epoch observation (see Table [TJ for each target was used. 
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number of free parameters is 4IV* + 4. We also only fit 
a circular aperture of radius 10 pixels around each star. 
From experimenting with different aperture sizes, we find 
this radius covered most of the star’s flux (the full width 
half maximum, FWHM, is about 5 pixels) while remain¬ 
ing small enough to avoid counting any remaining bad 
pixels in the background. We also explored some alter¬ 
native PSF fitting schemes using a smaller sample from 
our survey. We tried a Moffat function combined with an 
elliptical Gaussian model, a purely Gaussian model, and 
a sine 2 model. From examining the Bayesian Informa¬ 
tion Criteria, we find that the best model is the Moffat 
function with a radially symmetric (circular) Gaussian 
component. 

We find the best fit parameters using a maximum like¬ 
lihood estimation routine. These best fit parameters de¬ 
termine an analytic form for our PSF model. We com¬ 
pute the flux ratio by integrating the best-fit PSF model 
over the same circular aperture used in the PSF fitting 
for each star. We use the difference in the stellar posi¬ 
tion parameters to compute the horizontal and vertical 
separation as projected onto the NIRC2 array. We then 
adjust these separations to account for the well character¬ 
ized distortion and rotation of the NI RC2 array u sing the 
astrometric corrections presented in lYelda et alJ (i2010h . 

We compute the flux ratio and corrected one¬ 
dimensional separations for each frame. Our best esti¬ 
mate of these measurements for each observation is the 
median value from all of the frames. We calculate the 
standard error on the mean and use that as our mea¬ 
surement error. Using the best estimate for the one- 
dime nsional separations and the corrected NIRC2 plate 
scale (lYelda et al.ll2010H . we then compute the projected 
on-sky separation p and position angle PA for each de¬ 
tected companion. The NIRC2 astrometric corrections 
include uncertainties on the distortion, plate scale, and 
the orientation of the NIRC2 array. Therefore, we in¬ 
clude these uncertainties in our total error budget for 
our measured p and PA. 

We complete the above analysis to determine the pho¬ 
tometry for all detected companions in all bandpasses (J, 
H, K', K s ) and find the best fitting flux ratio between pri¬ 
mary and companion stars in each band, reported as a 
difference in magnitudes in Table [3] We then compute 
the apparent magnitudes of the companion stars in all 
bands as well as the colors of the companion stars. In 
order to do this, we require the apparent magnitudes of 
the prima ry star, w hich we o btained from the 2MASS 
catalog dSkrutskie et al.1 l2006fl . Table 2] shows the full 
multiband photometry of our detected companion stars. 
For our astrometric analysis, we only use K-band data 
and report the best fit separations and position angles in 
Tables [5] 

3.2. Common proper motion confirmation 

Our next step is to determine whether or not our can¬ 
didate companions share common proper motion with 
the primary star, indicating that they are bound com¬ 
panions rather than background sources in the same line 
of sight. If the detected companion is actually a very 
distant background star, it will remain effectively sta¬ 
tionary while the closer primary star moves across the 
sky as dictated by its parallax and proper motion. We 
would therefore expect that a background object would 


display a time-varying separation and position angle rel¬ 
ative to the primary star, while a bound companion will 
maintain a constant separation and position angle. 

We calculate the “background track” (i.e. the evolu¬ 
tion of the companion’s separation and position angle 
as a function of time if it was a background object) as 
shown in Figures [2] and [31 We compute the primary 
star’s parallactic motion using the celestial coordinates 
of the primary star from the SIMBAD database and the 
Earth ephemerides from the JPL Horizons service. The 
primary star’s proper motion and uncertainties are also 
taken from the SIMBAD database. When determining 
the background tracks, we account for uncertainties in 
the primary star’s celestial coordinates, proper motion, 
and parallax in addition to our measurement uncertain¬ 
ties in separation and position angle. We run a Monte 
Carlo routine to calculate the uncertainty in the back¬ 
ground tracks. The 68% and 95% confidence regions for 
the separation and position angle evolution are shown 
as shaded regions in Figures [2] and [31 We use the mea¬ 
surement with the smallest uncertainty in separation and 
position angle as the starting point for our track. 

After creating these background tracks, we next over¬ 
plot the measured companion separation and position 
angle at each epoch from Tables |1] [Gl and [7] Several 
of our candidate companions were detected in previous 
imaging surveys; when available, we also show the sep¬ 
arations and position angles from these earlier studies. 
We provide a complete list of these previously published 
detections in Tables [6] and H a nd d iscuss individual sys¬ 
tems in more detail in section 14.21 In Figures [2] and [31 
measurements from this study are plotted as circles while 
measurements from other studies are plotted as squares. 

Based on these measurements, we conclude that all 
but two of our detected candidate companions must be 
bound to the host star. The exceptions are the sec¬ 
ond candidate companion seen around HAT-P-7 and the 
candidate companion seen around HAT-P-15. Our fol¬ 
lowup observations determined these candidate compan¬ 
ions to be background objects. Due to the small pro¬ 
jected physical and angular separations of our compan¬ 
ion candidates, our result that the majority of our can¬ 
didates are physically bound companions is consistent 
with other direct imaging surveys th at confirm associa¬ 
tion via multi-epoch det ections (e.g. lEggenberger et all 
l2007t I Bowler e t all 1201 5fl or vi a galactic crowding esti¬ 
mates (eigTlAdams^Ol] I2012T) . We report confirmed 
common proper motion companions in Figure [2] and Ta¬ 
bles [5] and E The two background objects are reported 
separately in Figure [3] and Table [7] We discuss each sys¬ 
tem individually in the Section [4] below. 

3.3. Masses and separation 

For each confirmed companion-primary pair, we com¬ 
pute a flux ratio based on the stars’ physical parameters, 

f 0 °°I(\,T p ,log g p )R(\)r 2 p d\ 
r / 0 °° -^(A,T C , log g c )R(\)r% dA ’ 1 j 

where the subscripts p and c refer to the primary and 
companion stars, respectively, /(A, T, log g) is the star’s 
specific intensity, T is the star’s effective temperature, 
log g is a measure of the stellar surface gravity, R{\) 
is the response function of the NIRC2 filter, and r is 
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Figure 2. For each object, the two panels show the separation (top) and position angle (bottom) of a detected companion star relative to the 
primary target star. The lines show the track of a background object, as computed from our observation with the smallest uncertainty. The 
dark- and light-gray shaded regions indicate the 68% and 95% confidence region for this track. Filled symbols indicate measured positions 
of companions while open symbols indicate the position the detected source would have if it were a background source. Measurements 
from this work are plotted as circles and shown in Table [5] while measurements from previous studies are plotted as squares and shown in 
Table [6] Detected companions can be ruled out as background objects if either their separation or position angle are inconsistent with the 
background track. The two triple systems, HAT-P-8 and WASP-12, have each of their companion candidates plotted separately. Because 
the individual analysis for WASP-12B and WASP-12C cannot conclusively confirm or rule out common proper motion, we also plot the 
astrometric measurements for the center of mass for the combined light of both companions as WASP-12BC. All objects in this figure are 
determined to be common proper motion companions to their primary star. 


the stellar radius. In order to determine I(X,T,\ogg), 
we i nterpolate the grid ded PHOENIX synthetic spec¬ 
tra dHusser et al.l [2013t for solar metallicities and com¬ 
position ([Fe/H]= 0 and [a/H]= 0). Although our com¬ 
panion host stars have [Fe/H] between -0.21 and 0.25, we 
found that this assumption changes our estimated com¬ 
panion temperatures by less than 5 K. 

We use previously published measurements for the pri¬ 
mary star’s mass, radius, effective temperature and dis¬ 
tances as listed in Tableland fit for the companion star’s 
effective temperature, T c by matching the observed flux 


ratio with the computed flux ratio from Equation [2] In 
order to calculate the flux from a companion with a given 
effective tem perature, we use the z ero-age main sequence 
models from iBaraffe et al.l (119981) to match each effec¬ 
tive temperature with a corresponding radius and sur¬ 
face gravity, and then calculate the corresponding flux 
ratio using Equation [2] After determining the best-fit 
companion temperature, we calculate the corresponding 
uncertainty as the sum in quadrature of the uncertainties 
contributed by the flux ratio measurement error and the 
reported primary star temperature, surface gravity, and 
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Figure 3. These plots are the same as Figure(2] except this Figure 
shows the two candidate objects that were determined to be back¬ 
ground objects. The measurements are reported in Table [7] The 
top plot shows the astrometric measurements for the second can¬ 
didate companion to HAT-P-7, with circles showing our data and 
squares showing the four mea s ureme nts for the “western compan¬ 
ion” reported in INarita et all (12012 1. The bottom plot shows the 
astrometric measurements for the candidate companion to HAT- 
P-15. The data rules out a bound companion conclusion for both 
of these candidates. 


radius. We determine uncertainties for these parameters 
from uncertainties in our flux ratio measurement and pri¬ 
mary star’s temperature, radius, and logg. We do not 
include any uncertainties introduced through use of the 
stellar model or the PHOENIX model spectra. 

Finally, we convert our measured projected on-sky sep¬ 
arations to projected spatial separations using the stel¬ 
lar distance. Unfortunately most of our stars don’t have 
measured parallaxes; for these systems we generally used 
a spectroscopic distance estimate based on the derived 
stellar properties combined with the star’s apparent mag¬ 
nitude. 

We calculate an estimated temperature for each candi¬ 
date companion using either the J, H, or K band photom¬ 
etry, and find that these three independent temperature 
estimates are consistent with a late-type main sequence 
star in each filter. This indicates that all of the detected 
companions have infrared colors consistent with their in¬ 
ferred effective temperatures. In Table [8] we report the 
error-weighted averages of the companion stellar parame¬ 
ters from all three bands as well as temperature estimates 
using data from individual bands. We conclude that our 
detected companion stars have colors consistent with a 
late-type main sequence star in the same system rather 
than being a distant early-type star in the background. 

3.4. Contrast curves 
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Figure 4. The grey lines show K band contrast curve for all 50 
Friends of Hot Jupiters AO targets are plotted. For each target, 
the curve with the greatest contrast was chosen. All companion 
stars are masked out for the computation of these curves. The 
points show the magnitude difference and separation for all de¬ 
tected candidate companions at all epochs, where filled points in¬ 
dicate physically associated companions (with confirmed common 
proper motion) and open points indicate background objects. 

We calculate contrast curves for our target stars as fol¬ 
lows. First, we measure the FWHM of the central star’s 
PSF in the stacked and combined image, taking the aver¬ 
age of the FWHM in the x and y directions as our refer¬ 
ence value. We then create a box with dimensions equal 
to the FWHM and step it across the array, calculating 
the total flux from the pixels within the box at a given 
position. The ler contrast limit is then defined as the 
standard deviation of the total flux values for boxes lo¬ 
cated within an annulus with a width equal to twice the 
FWHM centered at the desired radial separation. We 
convert our absolute flux limits to differential magnitude 
units by taking the total flux in a box of the same size 
centered on the peak of the stellar point spread function 
and calculating the corresponding differential magnitude 
at each radial distance. We show the resulting 5cr average 
contrast curve for these observations in Figure [f] 

4. SYSTEMS WITH DETECTED COMPANIONS 
4.1. New physically bound companions 

Seven out of our 17 targets with detected companions 
(HAT-P-10, HAT-P-14, HAT-P-16, HAT-P-24, HAT-P- 
33, WASP-3, and WASP-14) have not been previously 
reported to have a directly imaged stellar companion. 
Here, we discuss each of these systems individually and 
report error-weighted averages of corresponding mea¬ 
surements in Tables [5] and [8] 

4.1.1. HAT-P-10 

HAT-P-10 hosts a transiting gas giant planet with 
a mass of 0. 5Mj„ n and an or bital semimajor axis of 
of 0.04 AU (jBakoset all l2009aih HAT-P-lOb does not 
have a published Rossiter-McLaughlin measurement con¬ 
straining its spin-orbit measurement, but its eccentricity 
is co nsistent with zero (jBakos et all l2009al iWest et all 
120091 ). We therefore placed this system in our control 
sample. We report a 0.36 Mq± 0.05 Mgstellar companion 
at a projected separation of 42 AU±2 AU. In the first pa¬ 
per o f the Friends of Hot Jupiter survey (jKnutson et all 
|20H . HAT-P-10 was found to have a detected radial 
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velocity trend that is consistent with the mass and sepa¬ 
ration of the directly imaged companion reported in this 
work. We show that HAT-P-10B is a common proper mo¬ 
tion companion. Although there are three other systems 
with both radial velocity trends and a directly imaged 
stellar companion, this is the only such target where the 
stellar companion might pl ausibly explain the observed 
trend (IKnutson et al.ll20l4lh 

4 . 1 . 2 . HAT-P-14 

iTorres et al.l (I2010T ) report a transiting gas giant planet 
around HAT-P-14 with a mass of 2.2 Mj up and an orbital 
semimajor axis of 0.06 AU. We placed this system in 
our misaligned sample because th is pl anet has an eccen¬ 
tricity of 0.12 ± 0.02 (IKnutson et al.l I2 014T) a nd a spin- 
orbit angle of —171° ± 5° (iWinn et al.ll2011l ). We find 
a 0.20 Mq± 0.04 Mq stellar companion to HAT-P-14 at 
a projected separation of 180 AU ± 10 AU. Our observa¬ 
tions in 2012 through 2014 confirm that HAT-P-14B is a 
common proper motion companion. 

4 . 1 . 3 . HAT-P-16 

HAT-P-16 hosts a transiting gas giant planet with 
a mass of 4.2M.T„ n and an or bit al semimajor a xis o f 
0.04 AU dBuchhave et al.l feOlOl h iMoutou et al.l (12011D 
measured a sp i n-orb it angle consistent with zero but 
IKnutson et alj (l2014f l measures a planet orbital eccen¬ 
tricity of 0.04 ± 0.01 so we placed this system in our 
misaligned sample. We report a 0.19 M 0 ± 0.03 M 0 stel¬ 
lar companion to HAT-P-16 at a separation of 160 AU ± 
10 AU. From our detections in February 2012 and Au¬ 
gust 2013, we are able to confirm that HAT-P-16B is a 
common proper motion companion. 

4 . 1 . 4 . HAT-P-24 

iKipping et all ( 2010 1 report a transiting gas giant 
planet around HAT-P-24 with a mass of 0.7Mj up and 
an orbital semimajor axis of 0.05 AU. They also find 
the planet’s eccentric ity to be consistent with zero and 
lAlbrecht et alj (I2012T) find evidence for a well-aligned or¬ 
bit, so we placed this target in our control sample. We 
report a 0.33 M 0 ± 0.03 M 0 stellar companion to HAT- 
P-24 at a separation of 2000 AU ± 100 AU. From our 
detections in February 2012 and May 2013, we are able 
to confirm that HAT-P-24B is a common proper motion 
companion. 


4 . 1 . 5 . HAT-P-33 

HAT-P-33 is a planetary system with a transiting gas 
giant planet with a m ass of 0.8 M.T 1]r ,and an o rbital semi¬ 
major axis of 0.05 AU (lHartman et al.ll2~011dh We placed 
this system in our control sample beca use H A T-P-33b 
has an eccentricity consistent with zero (lHartman et al.l 
l2011cl ) and no published spin-orbit angle measurement. 
We report a 0.55 M^O.OS M 0 stellar companion to HAT- 
P-33 at a separation of 119 AU ± 3 AU. From our detec¬ 
tions in February 2012 and March 2013, we are able to 
confirm that HAT-P-33B is a common proper motion 
companion. 


0.03 AU (jPollacco et al.l l2008lh WASP-3b has an ec¬ 
centricity consiste nt with zero (|Pollacco et al.1 120081 : 
iGibson et al.ll2008f) and a measured sp in-orbit angle con¬ 
sistent with zero dTrimthi et al.ll201 0l so we placed this 
target in our control sample. We report a 0.108 M 0 ± 
0.006 M 0 stellar companion to WASP-3 at a separation 
of 300 AU ± 20 AU. From our detections in June 2012, 
July 2012, and May 2013, we are able to confirm that 
WASP-3B is a common proper motion companion. 


4 . 1 . 7 . WASP-14 

Lloshi et all ( 2009 !) report a transiting gas giant planet 
around WASP-14 with a mass of 7.3Mj up and an or¬ 
bital semimajor axis of 0.04 AU. This system is in our 
misaligned sampl e because the p l anet h as an eccentricity 
of 0.082 ± 0.003 (IKnutson et al.l I2014T) a nd a spin-orbit 
angle of —33° ± 7° ( Johnson et al.ll2009lh We report a 
0.33 M 0 ± 0.04 M 0 stellar companion to WASP-14 at a 
separation of 300 AU ± 20 AU. From our detections in 
July 2012 and March 2013, we are able to confirm that 
WASP-14B is a common proper motion companion. 


4 . 2 . Companions reported in previous studies 

The remaining ten stellar companions were detected 
in previous studies, although not all systems were con¬ 
firmed to be physically bound. We discuss each system 
individually below, with measurements reported in Ta¬ 
ble El 


4 . 2 . 1 . HAT-P-7 


HAT-P-7 hosts a transiting gas giant planet with 
a mass of 1.8Mjup and an orbital semimajor axis of 
0.04 AU (lPa.1 et al.1 120081 b HAT-P-7b has a mea sured 
eccentricity consistent with zero ( EaLet_aP I2Q08 ) but 
its spin orbit ang le is 160° ± 40° ( Albrecht et al.l 20121 
IWinn et al.l 120091) so we placed this target in our mis¬ 
aligned sample. Our reported stellar companion to HAT- 


P-7 is consistent w ith the seven previous detections by 
iFaedi et ahl (1201 31): iNarita et al.1 (|2010bl b iBergfors et al.1 
(|2013l) : INarita et al.l (120121) b etween 2009 and 2012. The 
study bv INarita et al.1 ( 20121) determined this detected 
stellar companion HAT-P-7B is common proper mo¬ 
tion companion. Our own detections in 2012 through 
2014 are also con s istent with this conclusion. In addi¬ 
tion, INarita et al.l (2012D also reported a second directly 
imaged stellar companion (a “western companion”) in 
four of their images (spanning the same dates as above) 
with a separation of 3.1" ± 0.1" and a position angle of 
265.9° ± 0.4°. They were unable to confirm that this ob¬ 
ject is bound to the primary. We were able image this 
second companion on three out of the four nights we ob¬ 
served HAT-P-7. Our analysis (Figure [3]) concludes that 
this second companion is actually a background object 
and therefore not physically bound to the other two stars. 
This target is one of four systems with both a directly 
imaged stellar compani on and a long term r adial veloc¬ 
ity trend, as reported in lKnutson et al„ : ( 20141 b However, 
our resolved stellar companion, at a projected separation 
of 1240 AU±70 AU, is too distant to explain the observed 
radial velocity trend. 


4 . 1 . 6 . WAS P-3 

WASP-3 hosts a transiting gas giant planet with 
a mass of 1.8Mj up and an orbital semimajor axis of 


4 . 2 . 2 . HAT-P-8 

iMancini et al.l (120131 1 report a transiting gas giant 
planet around HAT-P-8 with a mass of 1.3Mj up and an 
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orbital semimajor axis of 0.04 AU. They also report the 
planet’s eccentricity to be consistent with zero and its 
spin- orbit angle has been measured to be consist ent with 
zero (iSimoson et al.lf20nl iMoutou et all 1201 ll) . Thus, 
we placed this system in our control sample. In our sur¬ 
vey, we find two stellar companions for HA T-P-8. These 
compa nions were previously detected by iBergfors et al.l 
(i2013l) in 2009 but their photometry did not resolve the 
two stars and they instead reported a single companion. 
We resolve these two stars in our Kec k images, which 
were p ublished in iBechter et al.l (12014D . iBechter et al.1 
(|2014li followed up on this target in 2012 and 2013 and 
find that both stellar companions has common proper 
motion with the primary. In this work we present a 
new set of observations, also taken in 2012 and 2013, 
which also indicate that HAT-P-8 is a physically bound 
triple system. We d o not use the unresol ved astrometric 
measurements from Bcrgfor s et all (120131) in our analysis 
becau se the resolved measurements from IBechter et al.1 
(|2014f) and this survey were enough to confirm common 
proper motion. 


4.2.3. HAT-P-30 

HAT-P-30 is a planetary system with a transiting gas 
giant planet with a mass of 0. 7M/inr>and an o rbital semi¬ 
major axis of 0.04 AU ([Johnson et al.l 120111 ). Although 
they measure the planet’s eccentricity to be consistent 
with zero, they also measure a spin-orbit misalignment 
of 74° ± 9°. We therefore placed this target in our 
misaligned sample. Our detected stellar companion is 
consis tent with the companion reported bv I Adams et al.l 
(120131) . observed between 2011 Oct 9-11. With only one 
detection, they were unable to confirm common proper 
motion. We found this companion in 2012 and followed 
up in 2014. Our data confirm that the com panion HAT - 
PdlOB is indeed bound. We do not use the lAdams et al.1 
(120131) measurement in this analysis as they did not re¬ 
port uncertainties on their measured separation and po¬ 
sition angle. 


4.2.4. HAT-P-32 

HAT-P-32 hosts a transiting gas giant planet with 
a mass of 0.9Mj lln and an orbi tal semimajor axis 
of 0.03 AU (lHartman et al.l 1201 Id) . HAT-P-32b has 
an ec centricity consistent with zero (jHartman et al.1 
1201: 3) but a mesaured sp in-orbit misalignment of 85° ± 
2° ( Albrecht et all 12012 ) so we placed this system in 
our misaligned sample. Our detected stellar compan- 
ion is consist ent with the stellar companion reported by 
lAdams et al.l (120131) . observed between 2011 Oct 9-11. 
As with HAT-P-30, they were unable to measure the 
proper motion for this companion. We found this stellar 
companion in 2012 and followed up in 2013. Our data 
confirm s that HAT -P- 32B is indeed bound. We do not 
use the lAdams et al.1 (120131) measurement in this analy¬ 
sis as they did not report uncertainties in the measured 
separation and position angle. This target is also one of 
four targets with both a directly imaged stellar compan¬ 
io n and a long term ra dial velocity trend, as reported 
in iKnutson et al.1 (|2014D . However, our resolved stellar 
companion, at a projected separation of 830 AU± 20 AU, 
is too distant to explain the observed radial velocity 
trend. 


4.2.5. TrES-2 


iBarclav et al.l ( 2012 ) report a transiting gas giant 
planet around TrES-2 with a mass of 1.4Mj up and an 
orbital semimajor axis of 0.04 AU. They a lso report a 
pl anet eccentricity consistent with zero and I Winn et al.l 
( 2008a ) report a spin-orbit angle consistent with zero so 
we placed this target in our control sample. Our detected 
stellar companion for TrES-2 is c onsistent with the thre e 
previous de t ection s repo rted bv IDaemgen et all (|2009l) . 
IFaedi et all (120131) . and IBergfors et ahl i 20 131) between 
2007 and 2009. These studies were not able to measure 
the proper motion of the companion in order to deter¬ 
mine whether or not it was bound to TrES-2. Our mea¬ 
surements from 2012 and 2013 show that TrES-2B is a 
common proper motion companion. 


4.2.6. TrES-4 

TrES-4 is a planetary system with a transiting gas gi¬ 
ant planet with a mass o f 0.9 M.T llrl and an o rbital semi¬ 
major axis of 0.05 AU dSozzetti et al.l 120091) . TrES-4b 
has a n ec centricity consistent with zero (|Sozzetti et al.l 
120091) and iNarita et al.1 (l2010al) measure a spin-orbit an¬ 
gle consistent with zero so we placed this system in 
our control sample. Our detected stellar companion for 
TrES-4 is c onsistent with the thre e previous de t ection s 
rep orted bv IDaemgen et all (|2009D . IFaedi et all (|2013l) . 
and IBergfors et ah (120131) between 2008 and 2009. Us¬ 
ing al l the data from the prior studies, Bcrgfors et al. 
dm! were able to confirm that TrES-4B is a common 
proper motion companion. Our measurements from 2012 
and 2013 support this assessment. 


4.2.7. WASP-1 

WASP-1 hosts a transiting gas giant planet with a mass 
of 0.8 Mj up (IKnutson et al.||2014|) and an orbital semima- 
i or ax is of 0.04 AU ( Torres et al.l 120081) . lAlbrecht, et all 
(120111) report that the eccentricity and spin-orbit angles 
are consistent with zero, so we placed this ta rget in our 
control sample. iCollier Cameron et al.l (120071) observed a 
stellar companion to WASP-1 with a separation of 4.7" 
north of the primary from observations made in 2006. 
They did not report any uncertainties nor did they pro¬ 
vide a numerical positional angle. Thus, we are unable to 
include this observation in our analysis. However, these 
values are consistent with the reported separation and 
position angle from our observations in 2012, 2013 and 
2014. Our measurements confirm that WASP-1B is a 
common proper motion companion. 


4.2.8. WASP-2 

There is a transiting gas giant planet around WASP- 
2 with a mass of 0.9Mj up (IKnutson et al.||2014[) and a n 
orbital semimajor axis of 0.06 AIJ (jTorres et al.ll2008D . 
Radial velocity and secondary eclipse measurements for 
this planet a re co nsistent with a circular orbit, but 
iTriaud et al.1 (120101) find that the planet has spin-orbit 

misalignment of 150°_ 2 qo . We therefore place this 
planet in the misaligned sample. A stellar compan¬ 
ion for WASP-2 was reported in six separate obser¬ 
vation s from 2006 to 2 011 b y ICollier Cameron et al.1 
(120071), IDaemgen et al.1 (120091) , IBergfors et al.l ( 2013D . 
and lAdams et alT l 20131) . These observations are all 






































































































































10 


Ngo et al. 


consistent with the companion we detect in our im¬ 
ages^ H owever, the m easurements reported in 2006 by 
iCollier Cameron et al.l (120071) did not include uncertain¬ 
ties nor a numerical value for the position angle, so we do 
not use this observation in our analysis. I n add ition, one 
observation from 2011 bv lAdarns et al.1 ( 2013 :) also did 
not include uncertaintie s so w e omit this measurement 
as well. iBergfors et al.1 (120131) combined their observa¬ 
tions with previous studies to conclude that the stellar 
companion must be bound to the primary. Our own mea¬ 
surements from 2012 and 2013 also indicate WASP-2B is 
a common proper motion companion. 

4.2.9. WASPS 

WASP-8 is a planetary system with a transiting gas gi¬ 
ant planet with a ma ss of 2.2 M.i lin and an orbital semima¬ 
j or ax is of 0.08 AU dQueloz et al.l 1201 0l) . Knutson et al.1 
(120141) reports the planet’s eccentricity to be 0.304±0.004 
and the discovery paper measures a spin-orbit misalign- 

■ o o 

ment of —123°_4 0 . Therefore, we placed this planet in 
the misaligned sample. This paper also reported a stel¬ 
lar companion for WASP-8 consistent with our own ob¬ 
servations in 2012 and 2013. They did not provide a 
date for their observation, so we do not use their mea¬ 
surement in our analysis. However, they note that this 
companion w as mentioned in the Washington Double 
Star Catalog (jMason et al.l l200ll) from 70 years ago at 
the same position. They conclude that WASP-8B is a 
common proper motion companion. Our data also rules 
ou t a background obj ect and supports the conclusion 
of iQueloz et al.l (l2010i) . This target is one of four tar¬ 
gets with a directly imaged companion and also a l ong 
term r adial velocity trend, as reported in iKnutson et al.l 
• 201 1:. However, our resolved stellar companion, at a 
projected separation of 380 AU ± 50 AU, is too distant to 
explain the observed radial velocity trend. 

4.2.10. WASP-12 

WASP-12 hosts a transiting gas giant planet with 
a mass of 1. 4M.T, ln and a n orbital semimajor axis of 
0.02 AU dHebb et al.1 120091). WASP -12b has eccentric¬ 
ity consistent with z ero (|Hebb et al.l 120091) . however, 
lAlbrecht et al.l (120121) measures a spin-orbit angle of 
59°i2o°- Thus, we placed this system in our mis¬ 
aligned sample. Stellar companion(s) for WASP-12 have 
been seen in 7 observat i ons fr o m 2009 to 2013 as re ¬ 
ported by Bergfors et al.1 (120131) ICrossfield et al.1 (|2012f ) . 
ISing et ahl ~( 20131) . and iBechter et al.1 (120 biff The first 
two studies were not able to resolve the two stellar com¬ 
panions for WASP-12 but did dete c t a s ingle source 
at the same position. IBergfors et al.1 (12 0131) noted that 
the source was elongated and ICrossfield et al.l (I2012D ob¬ 
tained a spectrum for the companion, which they used to 
estimated its effective temperature and surface gravity. 
They concluded that the observed source must be a fore¬ 
ground object as it was otherwise signi ficantly b righter 
than expected for its spectral type. ISing et all ( 1201 3) 
was the first to report images of the two resolved stellar 
companions but did not provide any astrometric mea¬ 
surements so we do not include this dat a point in our 
analys is. The two observations reported in lBechter et al.l 
(120141) are part of the Friends of Hot Jupiters survey. We 
independently analyzed the data and our results are con¬ 


sistent with the measurements reported in lBechter et al.l 
(120141) . Our measurement uncertainties are lower as we 
account for the NIRC2 distortion; thus we use the mea¬ 
surements reported in this work for our analysis. We 
also recently di scovered a more pre cise distance estimate 
for this s tar inlTriaud et ahl (120141) . which we use to re¬ 
place the IBergfors et al.l ( 20131) value from our previous 
papers. Despite these improvements, we find that our 
2012 and 2013 data points are insufficient to rule out the 
possibility of a background source (see astrometric plots 
for WASP-12B and WASP-12C separately in Figure [2] 
these plots do not include the combined light astrometric 
measurements). However, when considering previously 
reported combined light measurements for this system 
our analysis indicate the center of mass of WASP-12BC 
has common proper mo tion with the primar y, consistent 
with our conclusions in lBechter et al.l 


4.3. Candidates determined to be background objects 

As discussed in Section 14.2.11 we found that the 
second cand i date companion to HAT-P-7 reported in 
iNarita et al.l (|2012l) and followed up in our survey was a 
background object (see Figure [3] and Table [7]), so HAT- 
P-7 is simply a binary system. 

We also detect a candidate companion to HAT-P-15 
on 2 February 2012 with AK' = 5.77 ± 0.05, a separa¬ 
tion of 7.100"±0.002" and position angle of 58.6° ±0.1°. 
Due to the large separation of this object to HAT-P-15, 
the object was only visible in one of our dither posi¬ 
tions in our first epoch, resulting in relatively large as¬ 
trometric uncertainties. We followed up with a second 
epoch on 3 October 2014 and did not find an object in 
the same position. However, we did find an object with 
AK S = 8.2 ± 0.9, a separation of 6.387" ± 0.008" and po¬ 
sition angle of 58.0° ± 0.1°. Physical association of these 
two detections are ruled out by our astrometric analysis 
(see Figure [3] and Table [7]). In addition, the location of 
these detections are also inconsistent with both detec¬ 
tions being the same object. Thus, we conclude that we 
imaged two different background or foreground objects 
in the two epochs and exclude this target from our list 
of bound companions. 


4.4. Non-detections 

For the remaining 32 targets, we did not find any 
candidate companions within the 5.5" NIRC2 field of 
view. For two of these targets, previous studies have 
found directly i maged companions at larger separations. 
iMugrauer et al.l (120141 ) reports a common proper motion 
companion to the south of HAT-P-4 with a separation 
of 91.8". IBergfors et al.l (I2013D reports a candidate com¬ 
panion to XO-3 with a separation of 6.059" ± 0.047" and 
a position angle of 296.7° ± 0.3°. This companion is also 
very faint (A z' = 8.22 ± 0.23 and Az' = 8.57 ± 0.24) and 
IBergfors et al.1 (I2013D note that it is unlikely that the de¬ 
tected object is a bound companion. This companion 
would not be visible in our survey due to its faintness 
and our survey’s FOV, given our 3-point dither pattern. 

5. COMPANION FRACTION 

We find bound stellar companions around 17 out of 
the 50 targets observed, corresponding to an overall raw 
companion fraction of 34% ± 7%. We find that 9 out of 
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Figure 5. For each target, we plot the weighted average of com¬ 
panion masses and projected separation from all epochs reported 
in Table[H] Only one point is plotted for each of the triple systems. 
We plot targets with a non-zero spin-orbit angle measured using 
the Rossiter-McLaughlin (RM) effect as black diamonds, targets 
with spin-orbit angles consistent with zero as red circles, and tar¬ 
gets without any spin-orbit angle measurement as open squares. 
The two objects without a RM measurement are HAT-P-lOb and 
HAT-P-33b; both have orbital eccentricities consistent with zero 
so we placed them in the “control” sample. Companions that lie 
above and to the left of the solid, dashed, and dotted lines are 
sufficiently massive and close enough to excite Kozai-Lidov oscil¬ 
lations and overcome general relativity pericenter precession for 
giant planets at 1 AU, 2.5 AU and 5AU, respectively. 

27 stars with planets on misaligned or eccentric orbits 
have stellar companions, yielding a raw stellar compan¬ 
ion fraction of 33% ± 9%. We also find that 8 out of 
23 stars with planets on well-aligned or circular orbits 
have stellar companions, corresponding to a raw stellar 
companion fraction of 35%±9%. Figure [5] plots the phys¬ 
ical separation of detected companions versus their mass 
ratios. Companions to stars with measured spin-orbit 
misalignment are shown in black, companions to stars 
with measured spin-orbit angles consistent with zero are 
shown in red, and companions to stars with no spin- 
orbit angle measurement (HAT-P-10 and HAT-P-33) are 
shown as open squares. We find no evidence for a differ¬ 
ence in the typical mass ratios or projected separations 
between misaligned and well-aligned targets. 

Figure [5] also plots lines to show the minimum com¬ 
panion mass necessary to excite Kozai-Lidov oscillations 
at a timescale short enough to overcome pericenter pre¬ 
cession due to general relativity. For these representa¬ 
tive lines, we use a central stellar mass of 1 M@, a plan¬ 
etary mass of 1 Mj up , and assume a circular orbit for 
both planet and companion star. We plot these lim¬ 
its for initial planet semimajor axis distances of 1 AU, 
2.5 AU and 5AU (solid, dashed and dotted lines, re¬ 
spectively). These lines are computed by equating the 
timescales for the Kozai-Lidov oscillation and general 
relativity pericenter precession, usin g Equations 1 and 
23 from lFabrvckv fc Tremaine ! (2007). These expressions 
scale as 1 — e 2 , so they are not strongly affected by our 
assumption of circular orbits. All companions that lie 
above and to the left of these lines are sufficiently mas¬ 
sive and close enough to excite Kozai-Lidov oscillations 
on the planet. We find that the majority of our detected 
companions could potentially excite Jupiter-mass planets 


that form within 5 AU. 

5.1. Survey incompleteness correction 

In order to make a reliable estimate of the companion 
fraction for our sample, we must correct for our survey 
incompleteness. For each of our 50 targets, we deter¬ 
mine our survey’s sensitivity to companions with various 
mass ratios and orbital semimajor axes. We create a 
50x50 grid of linearly spaced bins in mass ratio and log¬ 
arithmically spaced bins in semimajor axis out to the 
distance that corresponds to a separation of 5.5". For 
each grid point, we generate 1000 fake companions with 
a mass and semimajor axis within the grid point lim¬ 
its. We draw the companion’s orbital eccentricity from 
a unifor m distribution derived f rom surveys of field star 
binaries (jRaghavan ct al.1 1201011 and we randomize the 
remaining orbital elements. Then, for each fake compan¬ 
ion, we compute the projected angular separation and 
brightness ratio and compare to that target’s 5er contrast 
curve to determine whether or not our survey would have 
been able to detect that fake companion. Thus, we com¬ 
pute our survey sensitivity for star i for companions in 
the grid point corresponding to mass ratio m r and semi¬ 
major axis a to be Di(m r ,a). 

We then compute the total average sensitivity, or sur¬ 
vey completeness, for star i as S t by taking the average of 
all Di (m r , a) values weighted by the fre quency of a com¬ 
panion with each m r and a: f(m r ,a). iRaghavan et alJ 
(] 2010 r > determines f(m ri a) to be a log-normal distribu¬ 
tion in period (which is a function of m r and a). We 
use this distribution to compute each target star’s sur¬ 
vey sensitivity as: 

_ / / Di(m r ,a)f{m r ,a)dm r d\na 
f f f(m r ,a)dm r dlna 

where we integrate over the range of our survey. We are 
sensitive to companions with periods between 10 4 days 
and 10 7 ' 5 days, which approximately corresponds to 
semimajor axes between 50 AU and 2000 AU. 

We use our estimate of survey completeness for each 
star to compute the likelihood L of obtaining our data 
set of N d detected companions out of N survey targets 
as 

N d N-N d 

l = n^ 7 ?) n w 

*=i i= i 

where 77 is the true companion fraction and the product 
sum over i is for the targets with a detected companion 
while the product sum over j is for the targets without 
a detected companion. We define companion fraction as 
the fraction of stars with at least one stellar companion 
within our survey range. 

We use the Affine-Invariant Markov Chain Monte 
Carlo scheme implemented by t he python package “em¬ 
cee” (j Goodman fc Wearel 1201(1 iForeman-Mackev et al.l 
12013! ) to determine the posterior probability distribution 
of 77 . Our prior on 77 is uniform between 77 = 0 and 77 = 1. 
In addition, we assume that we are 100% complete for 
targets where we have detected at least one companion 
because we expect all triple or higher order systems to 
be hierarchical. That is, we set all S t = 1.0 but still com¬ 
pute Sj as described in Equation [3j This assumption is 
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Figure 6. A comparison of the posterior probability density of 
the fraction of stars with at least one companion, 77 , for our total 
sample, our misaligned sample, our contro l sample, and for sola r 
type field stars in the solar neighborhood JRaghavan et aT 1 l 2010 l) . 
Only companions with periods between 10 4 days and 10 75 days are 
considered. 


supported by the two hierarchical syste ms detected by 
our s urvey and by previous studies (e.g. lEggleton et al.1 
l2007lb 

We find that the companion fraction is 7 ?t = 49% ± 9% 
for our total sample, r)M = 48% ± 12% for our mis¬ 
aligned sample, and rjc = 51% ± 13% for our control 
sample. These posterior distributions are shown in Fig¬ 
ure [6j These fractions are consistent with one another 
and we therefore conclude that there is no evidence for 
a correlation between the presence of a stellar compan¬ 
ion and the orbital properties of the transiting gas giant 
planet. 

We compare the companion fraction for hot (T e g > 
6200 K) and cool primary stars and find that hot stars 
have a higher companion rate, at the 2.9 a level. We find 
the companion fraction for hot stars to be 75%±14% and 
the companion fraction for cool stars to be 34% ± 10%. 
This difference is consistent with surveys of stellar mul¬ 
tiplicity, which indicate that more massive stars hav e 
higher binary fractions (e.g. iDuchene fc Krausl l2013h . 
Surveys of hot Jupiter obliquities indicate that mis¬ 
aligned hot Jupiters are preferentially located around 
stars hotter than 6200 K, which may be due to more ef¬ 
ficient tidal damping of primordial spin-or bit misalign¬ 
ment s in systems with cooler host stars (iSchlaufmanl 
120101 : iWinn et~aLl l2010t lAlbrecht et al .1 fJOTjfb We test 
this hypothesis by recalculating the companion fractions 
for our two sub-samples using only stars hotter than 
6200 K, where planets should still retain their primordial 
spin-orbit orientations. We find that hot stars in our mis¬ 
aligned sample have a companion fraction of 59% ± 17%, 
while hot stars in our control sample have a companion 
fraction of 83% ± 14%. These fractions are consistent at 
the 1.7 a level. 

We also consider the companion fractions for the 35 
stars in our sample with published measurements of spin- 
orbit alignment from the Rossiter-McLaughlin effect. We 
find the companion fraction for stars with a non-zero 
spin-orbit angle to be 73% ± 15% while the companion 
fraction for stars with a spin-orbit angle consistent with 


zero is 53% ± 14%. These fractions are consistent at the 
1.4 a level. 

When we examined these same targets using the ra¬ 
dial velocity technique we found that 51% ± 10% had 
companion s with masses betwee n 1 — 13 Mj up on long pe¬ 
riod orbits dKnutson et al.ll2014f) . Aside from HAT-P-10, 
where the directly imaged stellar companion is consistent 
with the measured RV acceleration, we find no evidence 
for any correlation between the presence of a stellar com¬ 
panion and a measured RV acceleration. Approximately 
1 /3 of our target stars have a detected stellar companion 
and 1 /3 have a radial velocity acceleration, so we would 
therefore expect to see approximately one-ninth of the 
stars in our sample with both types of companions. We 
detect four such systems including HAT-P-10. We there¬ 
fore conclude that the rates of stellar and radial velocity 
companions appear to be independent of one another and 
sum the two in quadrature to obtain a combined stellar 
and planetary occurrence rate of 72% ± 16%. 


5.2. Comparison with other direct imaging surveys 

We compare this companion fraction to the results 
from other surveys for stellar companions to planet hosts 
and field stars. These surveys are summarized in Ta¬ 
ble O Our survey results are consistent with previous di¬ 
rect imaging surveys fo r stellar companions t o transiting 
gas g i ant planet hosts (iDaemgen et al.lf2 009; Fa cdi et all 
I2013t IBergfors et al.l 120131 : lAdams et al.l I2013H . These 
previous surveys had s mall sample sizes; the l argest sam¬ 
ple size was 21 stars (|Bergfors et al.l 1201311 . Our re¬ 
sults are also consistent with the result of infrared di¬ 
rect imaging surveys of Kepler t ransit i ng planet can¬ 
didat e hosts (lAdams et ah 120121 120131: iDressing et al.l 
IMi ILillo-Box et al.l 120121 [Mi l Wang et al.l |2014bf) 7 
Except for lWang et al.l ~ 2014blf and lHorch et al.l (12014 ), 
none of the above surveys quantify their completeness 
and their numbers are closer to our uncorrected com¬ 
panion fraction. In addition, these surveys have quite 
widely varying sensitivities that may contribute to the 
large observed scatter in measured companion fraction. 

Figure Q compares the companions reported in this 
work with these diffraction-limited near-infrared AO sur¬ 
veys. As companion masses are not always reported, we 
plot the difference in magnitude versus angular separa¬ 
tion for all reported companions. Our survey is primarily 
sensitive to companions within 5" while the others con¬ 
sider companions at larger separations. In addition, our 
survey is able to detect companions that are approxi¬ 
mately 2 magnitudes fainter than these previous studies 
at separations less than 2". Because the Kepler candi¬ 
date host stars are on average significantly more distant 
than our target stars, these surveys were unable to dis¬ 
tinguish between background objects and bound stellar 
companions. This also means that the Kepler surveys 
are less sensitive to stellar companions at small projected 
physical separations. We also limit our survey to systems 
with short-period transiting gas giant planets, while the 
Kepler planet candidate sample is dominated by much 
smaller planets, many of which are in compact multi- 
planet systems. 

The direct imaging surveys conducted at visible 
wavelengths (lLaw cd all 12014 iGilliland et al.1 120151 : 
iHorch et al.l 120141) report smaller companion fractions 
than our survey and other NIR surveys. This is not 
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Figure 7. The difference in magnitude (J or K s ) for our con¬ 
firmed physically associated companions compared to companion 
candidates from other near-infrared direct imaging surveys and 
a representative 5<r contrast curve from our survey. Our survey 
is more sensitive to faint companions at small separations. The 
published surveys shown here target Kepler planet candidate host 
stars, which are on average three times more distant than our tar¬ 
gets. Thus, our survey finds companions at smaller projected phys¬ 
ical separations as well. 

surprising, as the majority of our detected companions 
are M stars that would be much fainter at optical wave¬ 
lengths. These surveys also have relatively small fields 
of view, which would also contribute to a lower detec¬ 
tion rate for companions as compared to the wide-held 
infrared surveys. 

We note that our results for the multiplicity rate of 
hot Jupiter host stars may also have implications for the 
ongoing debate about the origin of the discrepancy be¬ 
tween hot Jupiter o ccurrence rates from radial velocity 
and transit surveys. IWang et~ahl (l2014af) point out that 
hot Jupiter occurrence rates in Kepler surveys are ap¬ 
proximately two to three times smaller than hot Jupiter 
occurrence rates in Doppler planet surveys. Because the 
Kepler survey does not filt er mult i-stellar systems from 
its target list, IWang et al.l (l2Q14alj suggest that one po¬ 
tential cause for the discrepancy may be that these com¬ 
panion stars are suppressing the formation of all planets, 
including hot Jupiters, in the Kepler sample. IWang et al.l 
: 201 lb : also finds that planets are 1.7 ± 0.5 times less 
likely to form in a system with a stellar companion within 
1000 AU. However, we find that half of all hot Jupiters 
are found in stellar binaries, indicating that stellar mul¬ 
tiplicity does not inhibit the formation of these systems 
for cases where the stellar companion is on an orbit ap¬ 
proximately 50-2000 AU. 

5.3. Comparison with field stars 

Finally, we compare our companion fraction to the 
measured multiplicit y rate for Solar-type fi eld stars. The 
most recent survey (IRaghavan et al.l 12010f) reports that 
44% ± 2% of Solar-type stars within 25 pc are in multiple 
star systems , which is in good agreemen t with previous 
results from IDuciuennov & Mavon (|1991f) . This value is 
corrected for survey completeness. This fraction includes 
companions with separations too small to resolve with 
direct imaging or too wide to be found within our sur¬ 
vey’s held of view. Thus, we only consider the compan¬ 


ions from IRaghavan et alJ (120101) with periods between 
10 4 days and 10' 5 days, which corresponds to physical 
separations of approximately 50 AU an d 20 00 AU. I f we 
select only the subset of binaries from IRaghavan et alJ 
(|2010f) that fall within this period range, we find that the 
fraction of held stars with stellar companions is 24%±1%. 
This is 2.8cr smaller than our estimated companion fre¬ 
quency for systems with short period gas giant planets, 
suggesting that the presence of a stellar companion in¬ 
creases the likelihood that a gas giant planet will migrate 
inward to a relatively short period orbit. 

6. SUMMARY 

We found nineteen companions around seventeen tar¬ 
gets, including two triple systems, HAT-P-8 and WASP- 

I2j _which we previously reported in iBechter et alJ 

(|2014f) . We measure the proper motions for all detected 
companions and confirm that they are physically bound 
to the planet-hosting primary stars. We report seven 
new multiple star systems with transiting giant plan¬ 
ets and provide follow-up observations for ten previously 
reported candidate multi-star systems. Our follow-up 
observations allowed us to confirm the bound nature 
of three previously detected candidate companions and 
were in good agreement with the conclusions of previous 
studies on the bound nature of the other seven compan¬ 
ions. We also determined t he second candidate compan¬ 
ion to HAT-P-7 found by INarita et al.l (1201211 to be a 
background object. For all systems, we provide updated 
astrometric measurements as well as estimated masses 
and physical separations for the observed companions. 
We find that most of the detected companions are mas¬ 
sive enough to excite Kozai-Lidov oscillations on giant 
planets forming within 5 AU. 

Our companion fraction is consistent with previous 
NIR direct imaging surveys of stellar companions around 
transiting planet hosts. We correct for survey sensitiv¬ 
ity and find that close-in transiting gas giant planetary 
hosts are approximately twice as likely to have at least 
one stellar companion with 50 AU < a < 2000 AU as 
compared to field stars, although the significance of this 
difference is only 2.8 a. We find that the companion frac¬ 
tion for systems hosting a transiting gas giant planet on 
a misaligned or eccentric orbit is indistinguishable from 
the companion fraction for systems hosting a planet on 
a well-aligned and circular orbit. This is consistent with 
other results that suggest hot Jupiters ma y not primari ly 
migrate via the Kozai-Lidov mechanism (jDawson et alJ 
120151 : iNaoz et ahll2012t iPetrovichl 1 201 ol) . We also recal¬ 
culate the companion fractions for our two sub-samples 
using only stars hotter than 6200 K, as it has been sug¬ 
gested that tidal evolution might be able to remove 
primordial spin-orbit misalignments for planets orbiting 
cooler stars. We find that there is no evidence for a 
correlation between the presence of a stellar companion 
and spin-orbit misalignment of the transiting hot Jupiter, 
in good agreement with our results for the full sample. 
Finally, we calculate the companion frequency for the 
overall sample as a function of host star temperature 
and find that stars hotter than 6200 K have a companion 
rate that is approximately two times larger (2.9(7 sig¬ 
nificance) than their cool cou nterparts, consistent wit h 
surveys of stellar multiplicity (iDuchcne fc Krausll2013l) . 

We conclude that the data are consistent with two 
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possible scenarios and discuss them here. First, stel¬ 
lar companions may not play a role in the determi¬ 
nation of hot Jupiter spin-orbit misalignments. This 
would be consi stent with the planet-planet scatter¬ 
ing scenario (e.g. iChatteri ee et al.ll 2008t iNaerasawa e t, al l 
20081: IWu fc Lithwickl 120111 iBeauge fc Nesvornvl l2012t 
Lithwick fe Wull 2014f). Howev er, our radial velocity sur- 
vev in iKnutson et al.l (12014H found no correlation be¬ 
tween the presence of a long term radial velocity accel¬ 
eration and the spin-orbit alignment of the inner tran¬ 
siting planet. Other proposed orbital obliquity excita¬ 
tion mechanisms that do not require a companion star 
include misalignment of the natal disk’s angular mo¬ 
mentum vector with respe ct to the stellar spin axis due 
to chaotic star formation dBate et al.1 120101: iThies et al~1 
120111 iFielding et all I2014H , magnetic warping torque s 
due to the primary star’s magnetic field (lLai et al.ll2011h . 
and m odulation of stella r surfa ces by internal gravity 
waves dRogers et al.ll2012L 1 20 13li . Alternatively, primor¬ 
dial gravitationally bound stellar companions may have 
acted to perturb the protoplanetary disks out of align¬ 
ment with their host stars at the epoch of star and planet 
formation (e.g. iBatvginl l2 012|: |Batyg n fc Adamd 120131 

iCrida fc Batvginl 12014 IStorch et al.l 2014H . However, 

dynamical processing by cluster evolution would have 
removed or exchanged these compani ons, diminis hing 
their c urrent observational signatures (jMalmbor g et al.1 
l2007bl i. In other words, the companions we observe to¬ 
day may not be the ones responsible for the facilitation 
of hot Jupiter misalignments. 

In this scenario, a majority of planetary systems should 
form from disks with random alignments, regardless of 
whether or not there is currently a stellar companion 
present. This prediction can be tested by measuring 
the spin-orbit alignment of a large number of copla- 
nar, multi-planet systems such as those detected by the 
Kepler survey. If a significant fraction of these sys¬ 
tems are misaligned even when no stellar companion is 
present, it would provide strong support for the ubiquity 
of primordial disk misalignments, which could also ex- 
plain the observed popu lation of misaligned hot Jupitcrs. 
iMorton & Wind 1 2014 1 present the framework for such 
a test. 

Whatever the favored scenario may be, it must also be 
consistent with the increased frequency of misaligned hot 
Jupiters found around hot stars. If a stellar companion 
is not responsible for misalignment, then tidal evolution 
of the star-planet pair, which proceeds at an enhanced 
rate aroun d cooler stars, could give r i se to the observed 
trend (e.g. lWinn et al.lfJO 1 0Hbaill20 1 2t I Valsecchi fe Rasiol 
l2014f) . However, this mechan ism has been recently criti¬ 
cized bv lRogers fc Lid (120131b who argue that such a pro¬ 
cess requires an unphysical set of assumptions, and would 
generally lead to a misalignment distribution that is in¬ 
consistent with the observed one. Alternatively, magnet¬ 
ically facilitated disk-star couplin g may cause cooler star s 
to realign with their disks dSnalding fc Batvginll2014bh . 
signaling constancy with the second scenario where pro¬ 
toplanetary disks are misaligned by transient stellar com¬ 
panions at early times. 

The apparent enhancement in the companion fraction 
for our sample of transiting planets versus that of field 
stars is suggestive, and may indicate that these compan¬ 
ions play a role in the migration process. Intriguingly, 


lLaw et ahl (l2014f) found tentative evidence that in the 
Kepler sample, short period gas giant planets are more 
likely to have stellar companions than their more distant 
counterparts, which is also consistent with the idea that 
these companions play a role in planet migration. In ad¬ 
dition, when considering the companion fraction f ound 
in our radial velocity survey (IKnutson et al.l [2014h . we 
find that the overall rate of both planetary and stellar 
companions in systems with close-in transiting gas giant 
planets is 72% ± 16%, suggesting that these systems fre¬ 
quently have companions that may interact dynamically 
with the short-period planet. 

In the future, we plan to survey a larger sample of plan¬ 
ets detected using the radial velocity technique, which 
will span a much broader range of semimajor axes than 
any of the transiting planet surveys. A number of 
sources show long-term Doppler accelerations indicating 
the presence of outer companions, and Keck AO imaging 
has been demonstrated as a successful technique for iden¬ 
tifying faint stellar and substell ar companions b ased on 
the existence of such “trends” (iCrenn et al.ll2012l ). These 
data should provide a more definitive test of the poten¬ 
tial correlation between multiplicity and orbital semirna- 
jor axis of the inner planet. We also plan to expand our 
sample of hot Jupiter systems with AO imaging, in order 
to reduce the uncertainties in our estimate of the stellar 
companion rate for these systems. 
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Table 1 

Summary of NIRC2 AO Observations 


Target 

N c 

UT Obs. Date 

Filter 

Array 

Tint 

N m 

-^stack 

GJ-436 

0 

2012 Feb 02 

K' 

256 

5 


18 

HAT-P-2 

0 

2012 May 29 

K' 

512 

13.3 


9 

HAT-P-4 

0 

2012 Feb 02 

K' 

1024 

15 


9 

HAT-P-6 

0 

2012 Jul 27 

K' 

1024 

9 


30 

HAT-P-7 

1 

2012 Jul 27 

J 

1024 

9 

10 

10 



2012 Jul 27 

H 

1024 

9 

10 

10 



2012 Jul 27 

K' 

1024 

9 

10 

10 



2013 May 31 

Ks 

1024 

9 

8 

8 



2013 Jun 22 

H 

1024 

9 

12 

12 



2013 Jun 22 

K s 

1024 

9 

12 

12 



2014 Jul 12 

K s 

1024 

15 

15 

15 

HAT-P-8 

2 

2012 Jul 27 

J 

1024 

9 

10 

15 



2012 Jul 27 

H 

1024 

9 

13 

15 



2012 Jul 27 

K' 

1024 

9 

15 

15 



2013 Aug 19 

H 

512 

12.5 

6 

6 



2013 Aug 19 

K' 

1024 

9 

12 

12 



2013 Aug 19 

Ks 

1024 

9 

12 

12 

HAT-P-10 

1 

2012 Feb 02 

J 

1024 

9 

9 

9 



2012 Feb 02 

K' 

1024 

10 

9 

9 



2013 Aug 19 

H 

1024 

9 

6 

6 



2013 Aug 19 

K' 

1024 

9 

12 

12 



2013 Aug 19 

K s 

1024 

9 

12 

12 

HAT-P-11 

0 

2012 Jul 27 

K' 

512 

5 


18 



2013 Aug 19 

K c 

1024 

9 


12 

HAT-P-12 

0 

2012 Feb 02 

K' 

1024 

15 


9 

HAT-P-13 

0 

2012 Feb 02 

K' 

1024 

9 


9 

HAT-P-14 

1 

2012 Jun 05 

K' 

512 

5 

33 

33 



2013 Mar 26 

H 

1024 

9 

5 

6 



2013 Mar 26 

K s 

1024 

10 

6 

6 



2014 Jul 07 

K s 

1024 

20 

5 

5 

HAT-P-15 

0 

2012 Feb 02 

K' 

1024 

15 


9 



2014 Oct 03 

K s 

1024 

20 


12 

HAT-P-16 

1 

2012 Feb 02 

J 

1024 

10 

14 

15 



2012 Feb 02 

K' 

1024 

15 

15 

15 



2013 Aug 19 

H 

1024 

18 

6 

6 



2013 Aug 19 

K s 

1024 

9 

12 

12 

HAT-P-17 

0 

2012 May 07 

K' 

1024 

0.9 


18 

HAT-P-18 

0 

2012 May 29 

K' 

1024 

30 


9 

HAT-P-20 

0 

2012 Feb 02 

K' 

512 

2.5 


18 



2013 Nov 17 

K c 

1024 

10 


12 

HAT-P-22 

0 

2012 Feb 02 

K' 

512 

10 


18 

HAT-P-24 

1 

2012 Feb 02 

K' 

1024 

15 

6 

9 



2014 May 13 

J 

1024 

10 

12 

12 



2014 May 13 

H 

1024 

20 

6 

6 



2014 May 13 

K s 

1024 

10 

12 

12 

HAT-P-26 

0 

2012 Feb 02 

K' 

1024 

15 


9 

HAT-P-29 

0 

2012 Feb 02 

J 

1024 

10 


9 



2012 Feb 02 

K' 

1024 

15 


9 

HAT-P-30 

1 

2012 Feb 02 

K' 

1024 

10 

9 

27 



2014 May 13 

J 

1024 

10 

12 

12 



2014 May 13 

H 

512 

10 

12 

12 



2014 May 13 

K s 

1024 

10 

12 

12 

HAT-P-31 

0 

2012 May 07 

K' 

1024 

8 


17 

HAT-P-32 

1 

2012 Feb 02 

J 

1024 

10 

9 

9 



2012 Feb 02 

K' 

1024 

15 

15 

15 



2013 Mar 02 

H 

1024 

5 

9 

9 



2013 Mar 02 

K s 

1024 

15 

15 

15 

HAT-P-33 

1 

2012 Feb 02 

J 

1024 

15 

8 

9 



2012 Feb 02 

K' 

1024 

15 

9 

9 



2013 Mar 02 

H 

1024 

10 

12 

12 



2013 Mar 02 

K s 

1024 

10 

12 

12 

HAT-P-34 

0 

2012 Jul 27 

K' 

1024 

9 


30 

HD-149026 

0 

2012 Jun 05 

K' 

256 

1 


39 



2013 Mar 26 

K s 

256 

9 


12 



2013 Jul 04 

H 

128 

10 


9 



2013 Jul 04 

K s 

256 

10 


12 

TrES-2 

1 

2012 Jul 27 

J 

1024 

9 

15 

15 



2012 Jul 27 

H 

1024 

9 

15 

15 



2012 Jul 27 

K' 

1024 

9 

15 

15 



2013 May 31 

H 

1024 

9 

12 

12 



2013 May 31 

K s 

1024 

9 

6 

6 

TrES-3 

0 

2012 Jul 27 

K' 

1024 

9 


30 

TrES-4 

1 

2012 Jul 27 

J 

1024 

9 

15 

15 



2012 Jul 27 

H 

1024 

9 

18 

18 
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Table 1 — Continued 


Target 

N c 

UT Obs. Date 

Filter 

Array 

Tint 

iVfit 

-^stack 



2012 Jul 27 

K' 

1024 

9 

18 

18 



2013 Jul 04 

H 

1024 

9 

9 

9 



2013 Jul 04 

K s 

1024 

10 

12 

12 

WASP-1 

1 

2012 Jul 27 

K' 

1024 

20 

20 

20 



2013 Aug 19 

H 

1024 

9 

6 

6 



2013 Aug 19 

K s 

1024 

12.5 

12 

12 



2014 Jul 12 

K s 

1024 

15 

15 

15 



2014 Oct 03 

Ks 

1024 

15 

12 

12 

WASP-2 

1 

2012 Jul 27 

J 

1024 

10 

15 

15 



2012 Jul 27 

H 

1024 

9 

15 

15 



2012 Jul 27 

K' 

1024 

9 

30 

30 



2013 Jun 22 

H 

1024 

9 

11 

11 



2013 Jun 22 

Ks 

1024 

9 

12 

12 

WASP-3 

1 

2012 Jun 05 

K' 

1024 

5 

30 

30 



2012 Jul 27 

J 

1024 

9 

14 

15 



2012 Jul 27 

H 

1024 

9 

15 

15 



2012 Jul 27 

K' 

1024 

9 

30 

30 



2013 May 31 

H 

512 

5 

5 

6 



2013 May 31 

K s 

1024 

9 

12 

12 

WASP-4 

0 

2012 Jul 27 

K' 

1024 

20 


15 

WASP-7 

0 

2012 Aug 29 

K' 

1024 

9 


48 

WASP-8 

1 

2012 Jul 27 

J 

1024 

9 

10 

15 



2012 Jul 27 

H 

1024 

9 

10 

15 



2012 Jul 27 

K' 

1024 

9 

20 

30 



2013 Aug 19 

K s 

512 

12.5 

8 

8 

WASP-10 

0 

2012 Jul 04 

J 

1024 

20 


9 



2012 Jul 04 

K' 

1024 

20 


9 

WASP-12 

2 

2012 Feb 02 

J 

1024 

15 

9 

9 



2012 Feb 02 

K' 

1024 

15 

9 

9 



2013 Mar 02 

K s 

1024 

10 

15 

15 

WASP-14 

1 

2012 Jul 27 

J 

1024 

9 

15 

15 



2012 Jul 27 

H 

1024 

9 

15 

15 



2012 Jul 27 

K' 

1024 

9 

33 

33 



2013 Mar 26 

K' 

1024 

9 

12 

12 



2013 Mar 26 

K s 

1024 

9 

12 

12 

WASP-15 

0 

2012 Feb 02 

K' 

1024 

15 


9 

WASP-16 

0 

2012 Feb 02 

K' 

1024 

15 


18 



2012 Jul 27 

H 

1024 

9 


15 



2012 Jul 27 

K' 

1024 

9 


33 

WASP-17 

0 

2012 May 07 

K' 

1024 

10 


18 

WASP-18 

0 

2012 Aug 28 

K' 

1024 

9 


32 

WASP-22 

0 

2012 Aug 26 

K' 

1024 

10 


9 

WASP-24 

0 

2012 Feb 02 

K' 

1024 

15 


9 



2012 Jun 03 

K' 

1024 

15 


12 

WASP-34 

0 

2012 Feb 02 

K' 

1024 

10 


18 

WASP-38 

0 

2012 Jun 04 

K' 

1024 

30 


4 

XO-2 

0 

2012 Feb 02 

K' 

1024 

10 


27 

XO-3 

0 

2012 Feb 02 

J 

512 

2.5 


9 



2012 Feb 02 

K' 

512 

2.5 


9 



2013 Mar 02 

H 

1024 

15 


12 



2013 Mar 02 

K' 

1024 

15 


9 

XO-4 

0 

2012 Feb 02 

K' 

1024 

9 


9 

XO-5 

0 

2012 Feb 02 

K' 

1024 

15 


16 


Note. — The N c columns indicates the number of confirmed common proper motion companions found for each target. The “Array” column 
denotes the horizontal width, in pixels, of the section of the NIRC2 detector used to capture the image. The array dimensions used in this survey 
were 1024x1024 (the full NIRC2 array), 512x512, 256x264, or 128x152. These dimensions are constrained by NIRC2’s readout software. The Tint 
column indicates the total integration time of a single exposure, in seconds. For targets with detected companions, the TVfit column indicates the 
number of exposures used to make our photometric and/or astrometric measurements of companions. The iV stac k column indicates the number of 
exposures used to compute the final stacked image for contrast curve measurements. The last column is only given for images taken in the K 7 , K s , 
or K c bandpasses because we only compute contrast curves in these bands. iVft t and iV stac k may differ because the companion may not be visible 
at all dither positions and/or poor conditions prevented acquisition of useful data. 
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Table 2 

Target stellar parameters 


Target 

N c 

Sample 

r e ff 

(K) 

M 

(M 0 ) 

log g 
(cgs) 

D 

(pc) 

References for... 

T M, log g D 

GJ 436 

0 

Misaligned 

3416 ± 54 

0.507 ± 0.070 

4.83 ±0.03 

10.14 ±0.24 

i 

i 

2 

HAT-P-2 

0 

Misaligned 

6290 ± 60 

1.36 ±0.04 

4.138 ±0.035 

125.3 ±13.1 

3 

3 

2 

HAT-P-4 

0 

Control 

5890 ± 67 

1.26 ±0.10 

4.12 ±0.03 

293.5 ± 19.4 

4 

5,6 2 

HAT-P-6 

0 

Misaligned 

6687 ± 71 

1.29 ±0.06 

4.188 ±0.035 

277.8 ± 19.1 

4 

7 

2 

HAT-P-7 

1 

Misaligned 

6259 ± 32 

1.361 ± 0.021 

4.98 ±0.13 

320+®° 

8 

8 

2 

HAT-P-8 

2 

Control 

6223 ± 67 

1.192 ± 0.075 

4.177 ±0.022 

230 ± 15 

4 

9 

2 

HAT-P-10 

1 

Control 

4974 ± 65 

0.83 ± 0.03 

4.56 ±0.02 

121.7 ±4.2 

4 

10 

10 

HAT-P-11 

0 

Misaligned 

4792 ± 69 

0.81 ± 0.03 

4.59 ±0.03 

38.0 ±1.3 

4 

11 

11 

HAT-P-12 

0 

Control 

4650 ± 60 

0.733 ± 0.018 

4.61 ±0.01 

139.1 ± 9.6 

12 

12 

2 

HAT-P-13 

0 

Misaligned 

5720 ± 69 

1.320 ±0.062 

4.070 ±0.020 

214 ± 12 

4 

13 

14 

HAT-P-14 

1 

Misaligned 

6671 ± 75 

1.418 ±0.054 

4.187 ±0.025 

205 ± 11 

4 

7 

15 

HAT-P-15 

0 

Misaligned 

5640 ± 69 

1.013 ± 0.043 

4.38 ±0.03 

190 ± 8 

4 

16 

16 

HAT-P-16 

1 

Misaligned 

6140 ± 72 

1.218 ±0.039 

4.34 ±0.03 

235 ± 10 

4 

17 

17 

HAT-P-17 

0 

Misaligned 

5345 ± 70 

0.857 ± 0.039 

4.52 ±0.02 

90 ±3 

4 

18 

18 

HAT-P-18 

0 

Control 

4790 ± 72 

0.770 ± 0.031 

4.57 ±0.04 

166 ±9 

4 

19 

19 

HAT-P-20 

0 

Misaligned 

4619 ± 72 

0.756 ± 0.028 

4.63 ± 0.02 

70 ±3 

4 

20 

21 

HAT-P-22 

0 

Control 

5367 ± 67 

0.916 ±0.035 

4.36 ± 0.04 

82 ±3 

4 

20 

21 

HAT-P-24 

1 

Control 

6329 ± 67 

1.191 ± 0.042 

4.27 ±0.04 

396 ± 20 

4 

22 

22 

HAT-P-26 

0 

Control 

5142 ±69 

0.816 ±0.033 

4.56 ±0.06 

134+1® 

4 

23 

23 

HAT-P-29 

0 

Control 

6086 ± 69 

1.207 ± 0.046 

4.34 ±0.06 

322+®® 

4 

24 

24 

HAT-P-30 

1 

Misaligned 

6304 ± 88 

1.242 ± 0.041 

4.36 ±0.03 

193 ±8 

25 

25 

25 

HAT-P-31 

0 

Misaligned 

6065 ± 100 

1.218 ±0.070 

4 26" 1 " 0 ' 11 
^• zo -0.13 

354 tit 

26 

26 

26 

HAT-P-32 

1 

Misaligned 

6207 ± 88 

1.160 ±0.041 

4.33 ±0.01 

283 ± 5 

27 

27 

27 

HAT-P-33 

1 

Control 

6446 ± 88 

1.375 ± 0.040 

4.15 ±0.01 

387 ±9 

27 

27 

27 

HAT-P-34 

0 

Misaligned 

6442 ± 88 

1.392 ± 0.047 

4.21 ±0.06 

257+^ 

28 

28 

28 

HD 149026 

0 

Control 

6103 ± 66 

1.345 ± 0.020 

4.189 ±0.020 

80.8 ±4.0 

4 

29 

2 

TrES-2 

1 

Control 

5850 ± 50 

0.94 ±0.05 

4.45 ±0.01 

220 ± 10 a 

30 

31 

2 

TrES-3 

0 

Misaligned 

5514 ± 69 

0.928 ± 0.038 

4.57 ±0.01 

258.5 ± 16.l a 

4 

32 

2 

TrES-4 

1 

Control 

6200 ± 75 

1.339 ± 0.086 

4.030 ± 0.033 

476 ± 26 a 

32 

7 

2 

WASP-1 

1 

Control 

6160 ± 64 

1.265 ± 0.054 

4.209 ±0.051 

380 ± 38 a 

4 

7 

2 

WASP-2 

1 

Misaligned 

5255 ± 71 

0.851 ± 0.050 

4.537 ±0.017 

140 ± 10 a 

4 

7 

2 

WASP-3 

1 

Control 

6375 ± 63 

1.20 ±0.01 

4.33 ±0.03 

220 ± 20 a 

4 

33 

2 

WASP-4 

0 

Control 

5540 ± 55 

0.927 ±0.056 

4.487 ±0.010 

280.9 ±31.1 

34 

7 

2 

WASP-7 

0 

Misaligned 

6520 ± 70 

1.317 ±0.072 

4.218 ±0.048 

140 ± 15 

34 

7 

35 

WASP-8 

1 

Misaligned 

5570 ± 85 

1.04 ±0.08 

4.41 ±0.09 

87 ±7 

34 

36 

2 

WASP-10 

0 

Misaligned 

4735 ± 69 

0.703 ± 0.070 

4 5i+ 0 - 06 
—0.05 

90 ±20 

4 

37 

37 

WASP-12 

2 

Misaligned 

6118 ±64 

1.38 ± 0.19 

4.159 ±0.024 

427 ± 90 a 

4 

7 

2 

WASP-14 

1 

Misaligned 

6462 ± 75 

1.350 ±0.121 

4.126 ±0.042 

160 ± 20 a 

4 

7 

2 

WASP-15 

0 

Misaligned 

6405 ± 80 

1.305 ± 0.051 

4.189 ±0.021 

256 ± 32 b 

36 

38 

See note b 

WASP-16 

0 

Control 

5630 ± 70 

0.980 ±0.054 

4.357 ±0.022 

174 ± 14 b 

36 

38 

See note b 

WASP-17 

0 

Misaligned 

6550 ± 100 

1.286 ± 0.079 

4.149 ±0.014 

476 ± 36 a 

39 

39 

2 

WASP-18 

0 

Control 

6368 ± 66 

1.274 ±0.060 

4.365 ±0.022 

122.6 ±6.7 a 

4 

7 

2 

WASP-22 

0 

Control 

5958 ± 98 

1.109 ± 0.026 

4.31 ±0.03 

300 ± 30 c 

40 

40 

41 c 

WASP-24 

0 

Control 

6107 ± 77 

1.184 ±0.027 

4.263 ±0.017 

332.5 ±23.8 a 

4 

42 

2 

WASP-34 

0 

Control 

5700 ± 100 

1.01 ± 0.07 

4.5 ±0.1 

120 ± 15 

43 

43 

43 

WASP-38 

0 

Misaligned 

6187 ± 77 

1.23 ± 0.04 

4.267 ±0.030 

110 ±20 

4 

44 

45 

XO-2 

0 

Control 

5377 ± 79 

0.924 ±0.173 

4.436 ± 0.042 

156.0 ±8.8 

4 

7 

2 

XO-3 

0 

Misaligned 

6759 ± 79 

1.213 ±0.066 

4.244 ± 0.041 

185.7 ±11.8 

4 

46 

2 

XO-4 

0 

Misaligned 

6297 ± 72 

1.32 ± 0.02 

4.18 ±0.07 

308.2 ± 19.6 

4 

47 

2 

XO-5 

0 

Control 

5370 ± 70 

0.88 ± 0.03 

4.31 ±0.03 

260 ± 12 

48 

48 

48 











References. — ( 

L) Ivon Braun et al. (120121); (2) Triaud et al. ( 20141); (3) Pal et al.l (20101): (4) Torres et al. (20121): (5)lKovacs et a .1 (120071): 

(6)IWinn et al.l (2011 

1) ; (7 ) Southworth ( 

20121); (8) IVan Even et al. (2012 1); (9) Mancini et al. (20131); (10)|Bakos et al.l d2009a|); (ll)IBakos et al. 

(120101); (12) 1 Hartman et al.| (120091); (13) 

Southworth et al. (2012al); (14) Bakos et al.l (2009bD; (15) Torres 

et al.l (120101); (16)|Kovacs et al.l d2010l); 

(17) Buchhave et al. 

(201 

0|); (18) IHoward et al. (l2012al); (19) Hartman et al.| ( 201 lal) ; (20) Bakos et a 

.1 (120111); (21) IBakos et al. (120101); (22) 

Kipping et al.| (12010 

; (23) Hartman et al.l (12011 b|) ; (24) Buchhave et al.| (20111); (25) IJohnson et al.| 

20 111); 

(26) 

Kipping et al. (20111); (27) 

Hartman et al.| (2011 3) ; (28)[Bakos et al.ld2012f); (29)ICarter et al.l (20091); (30)ISozzetti el 

al.1 120071): (31') 

Barclav et a 

. (120121); (32)|Sozzetti et al. 

(20091): (33) IMiller et al. 

120101); (34) Maxted et al.| (20111); (35) IHellier et al.l (120091); (36) IDovle et al. (20131); 

(37) IChristian et al. (20091); (38) 

Southworth et al.l d2013l); 

(39) ISouthworth et al. (!2012b|); (40) [Anderson et al. (20111); (41) Maxted et al.l (120101); (42) IStreet et al.l (120101); (43) 

Smallev et alj d2011|); (44) Brown et al. (201 2ft ; ( 45) IB 

axros et al.| d2010l); (46) IWinn et al.l (|2008b|); (47) McCullough 

et al. (2008D; (48) Pal et al. 


(120091) 


Note. — N c is the number of confirmed common proper motion companions found for each target. “Sample” shows whether we placed the 
target in the “misaligned” (planets on misaligned or eccentric orbits) or the “control” (planets on well-aligned and circular orbits) subsamples. 
With the exception of GJ 436 and HAT-P-2, our stars do not have directly measured parallax estimates. We take our distance estimates from 
the referenced papers, which combine estimates of the stellar surface gravity, effective temperature, and metallicity from high-resolution optical 
spectroscopy and (in some cases) constraints on the stellar density from fits to the transit light curve in order to constrain the mass, radius, and 
age of the host star. The measured apparent magnitudes in V , J, H , and K bands can then be used to estimate the distance to the star given 
these known properties. The final three columns lists references for our temperature; mass and log g ; and distance measurements, respectively. 

a This distance estimate has been updated from the values used in the first “Friends of Hot Jupiters” paper jKnutson et al1l2014D . The new 
distances are consistent with the previously used values except for WASP-12. 

k We could not find distance measurements for these targets so we estimated the distance based on relative brightness of the target to a reference 
star and used the same relative error. For WASP-15, we used WASP-14 as a reference and for WASP-16, we used WASP-8 as a reference. 

C A distance was reported with no uncertainty, so we use a conservative estimate of 10%, based on uncertainties of other stars. 
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Table 3 

Flux ratio measurements of confirmed common proper motion companions 


Companion 

UT Obs. Date 

AJ 

AH 

AK' 

AK 

: s 

HAT-P-7B 

2012 

Jul 27 

6.562 ± 0.075 

5.887 ± 0.078 

5.920 ±0.058 



HAT-P-7B 

2013 

May 31 




5.948 ± 

0.040 

HAT-P-7B 

2013 

Jun 22 


5.926 ± 0.093 


6.04 ± 

0.12 

HAT-P-7B 

2014 

Jul 12 




5.827 ± 

0.058 

HAT-P-8B 

2012 

Jul 27 

6.59 ± 0.12 

5.671 ±0.054 

5.874 ± 0.063 



HAT-P-8C 

2012 

Jul 27 

7.16 ± 0.15 

6.352 ± 0.065 

6.447 ± 0.057 



HAT-P-8B 

2013 

Aug 19 


5.597 ± 0.036 

5.841 ±0.073 



HAT-P-8C 

2013 

Aug 19 


6.147 ±0.077 

6.310 ±0.057 



HAT-P-8B 

2013 

Aug 19 


5.597 ±0.036 


5.726 ± 

0.075 

HAT-P-8C 

2013 

Aug 19 


6.147 ±0.077 


6.201 ± 

0.094 

HAT-P-10B 

2012 

Feb 02 

2.656 ± 0.017 


2.960 ±0.022 



HAT-P-10B 

2013 

Aug 19 


2.448 ± 0.014 

2.763 ±0.034 



HAT-P-10B 

2013 

Aug 19 


2.448 ± 0.014 


2.733 ± 

0.045 

HAT-P-14B 

2012 

Jun 05 



5.633 ±0.033 



HAT-P-14B 

2013 

Mar 26 


5.237 ±0.086 


5.647 ± 

0.096 

HAT-P-14B 

2014 

Jul 07 




5.844 ± 

0.034 

HAT-P-16B 

2012 

Feb 02 

5.421 ± 0.039 


5.530 ±0.021 



HAT-P-16B 

2013 

Aug 19 


5.147 ±0.047 


5.382 ± 

0.061 

HAT-P-24B 

2012 

Feb 02 



4.45 ±0.12 



HAT-P-24B 

2014 

May 13 

4.658 ± 0.080 

4.163 ±0.052 


3.811 ± 

0.010 

HAT-P-30B 

2012 

Feb 02 



3.134 ±0.054 



HAT-P-30B 

2014 

May 13 

3.4304 ± 0.0063 

3.1153 ±0.0053 


2.975 ± 

0.041 

HAT-P-32B 

2012 

Feb 02 

4.148 ± 0.028 


3.733 ±0.042 



HAT-P-32B 

2013 

Mar 02 


3.668 ± 0.071 


3.500 ± 

0.073 

HAT-P-33B 

2012 

Feb 02 

4.137 ±0.054 


3.938 ±0.085 



HAT-P-33B 

2013 

Mar 02 


3.469 ± 0.033 


3.415 ± 

0.057 

TrES-2B 

2012 

Jul 27 

2.896 ± 0.020 

2.4515 ± 0.0073 

2.422 ±0.019 



TrES-2B 

2013 

May 31 


2.4628 ± 0.0085 


2.523 ± 

0.036 

TrES-4B 

2012 

Jul 27 

3.528 ± 0.012 

3.177 ±0.013 

3.301 ±0.046 



TrES-4B 

2013 

Jul 04 


3.1841 ±0.0084 


3.447 ± 

0.024 

WASP-1B 

2012 

Jul 27 



4.551 ±0.077 



WASP-1B 

2013 

Aug 19 


4.766 ± 0.046 


4.840 ± 

0.055 

WASP-1B 

2014 

Jul 12 




4.858 ± 

0.064 

WASP-1B 

2014 

Oct 03 




4.561 ± 

0.100 

WASP-2B 

2012 

Jul 27 

2.880 ±0.010 

2.7343 ± 0.0061 

2.871 ±0.028 



WASP-2B 

2013 

Jun 22 


2.7254 ± 0.0075 


2.687 ± 

0.017 

WASP-3B 

2012 

Jun 05 



6.528 ±0.049 



WASP-3B 

2012 

Jul 27 

7.27 ±0.12 

6.683 ± 0.079 

6.641 ±0.052 



WASP-3B 

2013 

May 31 


6.550 ±0.091 


6.552 ± 

0.027 

WASP-8B 

2012 

Jul 27 

4.492 ± 0.062 

3.134 ±0.056 

2.550 ±0.021 



WASP-8B 

2013 

Aug 19 




2.560 ± 

0.018 

WASP-12B 

2012 

Feb 02 

3.711 ± 0.076 


3.32 ±0.11 



WASP-12C 

2012 

Feb 02 

3.672 ± 0.021 


3.577 ± 0.036 



WASP-12B 

2013 

Mar 02 




3.286 ± 

0.030 

WASP-12C 

2013 

Mar 02 




3.182 ± 

0.017 

WASP-14B 

2012 

Jul 27 

5.207 ±0.020 

4.788 ± 0.024 

4.984 ± 0.025 



WASP-14B 

2013 

Mar 26 



4.765 ±0.057 



WASP-14B 

2013 

Mar 26 




4.701 ± 

0.054 


Note. 


AX is the difference in magnitude between the companion and primary stars in the X filter. 
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Table 4 

Multi-band photometry of confirmed common proper motion companions 


Companion 

UT Obs. Date 

K 

rrij 

m H 

m K 

J-K 

H - K 

J-H 

HAT-P-7B 

2012 

Jul 27 

K' 

16.117 ±0.075 

15.231 ±0.078 

15.254 ± 0.058 

0.863 ± 0.095 

-0.023 ±0.097 

0.89 ±0.11 

HAT-P-7B 

2013 

May 31 

K s 



15.282 ± 0.040 




HAT-P-7B 

2013 

Jun 22 

K s 


15.270 ±0.093 

15.37 ± 0.12 


-0.10 ±0.15 


HAT-P-7B 

2014 

Jul 12 

K s 



15.161 ±0.058 




HAT-P-8B 

2012 

Jul 27 

K' 

15.81 ± 0.12 

14.675 ± 0.054 

14.827 ±0.063 

0.98 ±0.14 

-0.152 ±0.083 

1.13 ±0.13 

HAT-P-8C 

2012 

Jul 27 

K' 

16.37 ± 0.15 

15.356 ±0.065 

15.400 ± 0.057 

0.97 ±0.16 

-0.044 ±0.087 

1.02 ±0.16 

HAT-P-8B 

2013 

Aug 19 

K' 


14.601 ± 0.036 

14.794 ± 0.073 


-0.193 ±0.082 


HAT-P-8C 

2013 

Aug 19 

K' 


15.151 ±0.077 

15.263 ± 0.057 


-0.113 ±0.096 


HAT-P-8B 

2013 

Aug 19 

K s 


14.601 ±0.036 

14.679 ± 0.075 


-0.078 ±0.083 


HAT-P-8C 

2013 

Aug 19 

K s 


15.151 ±0.077 

15.154 ±0.094 


-0.00 ±0.12 


HAT-P-10B 

2012 

Feb 02 

K' 

12.671 ± 0.017 


12.381 ± 0.022 

0.290 ± 0.027 



HAT-P-10B 

2013 

Aug 19 

K' 


12.008 ±0.014 

12.184 ± 0.034 


-0.176 ±0.037 


HAT-P-10B 

2013 

Aug 19 

K s 


12.008 ±0.014 

12.154 ±0.045 


-0.146 ±0.047 


HAT-P-14B 

2012 

Jun 05 

K' 



14.484 ± 0.033 




HAT-P-14B 

2013 

Mar 26 

K s 


14.164 ±0.086 

14.498 ± 0.096 


-0.33 ±0.13 


HAT-P-14B 

2014 

Jul 07 

K s 



14.695 ± 0.034 




HAT-P-16B 

2012 

Feb 02 

K' 

15.271 ± 0.039 


15.083 ± 0.021 

0.189 ±0.045 



HAT-P-16B 

2013 

Aug 19 

K s 


14.770 ± 0.047 

14.935 ±0.061 


-0.165 ±0.077 


HAT-P-24B 

2012 

Feb 02 

K' 



14.99 ± 0.12 




HAT-P-24B 

2014 

May 13 

K s 

15.455 ± 0.080 

14.752 ± 0.052 

14.354 ± 0.010 

1.101 ±0.081 

0.398 ± 0.053 

0.703 ± 0.096 

HAT-P-30B 

2012 

Feb 02 

K' 



12.285 ± 0.054 




HAT-P-30B 

2014 

May 13 

K s 

12.8724 ± 0.0063 

12.3353 ± 0.0053 

12.126 ±0.041 

0.746 ±0.042 

0.209 ± 0.042 

0.5371 ±0.0082 

HAT-P-32B 

2012 

Feb 02 

K' 

14.399 ± 0.028 


13.723 ± 0.042 

0.676 ± 0.050 



HAT-P-32B 

2013 

Mar 02 

K s 


13.692 ±0.071 

13.490 ± 0.073 


0.20 ±0.10 


HAT-P-33B 

2012 

Feb 02 

K' 

14.400 ± 0.054 


13.942 ± 0.085 

0.46 ±0.10 



HAT-P-33B 

2013 

Mar 02 

K s 


13.530 ±0.033 

13.419 ± 0.057 


0.111 ±0.066 


TrES-2B 

2012 

Jul 27 

K' 

13.128 ±0.020 

12.3715 ±0.0073 

12.268 ± 0.019 

0.859 ± 0.028 

0.103 ±0.021 

0.756 ±0.021 

TrES-2B 

2013 

May 31 

K s 


12.3828 ± 0.0085 

12.369 ± 0.036 


0.014 ±0.037 


TrES-4B 

2012 

Jul 27 

K' 

14.111 ± 0.012 

13.527 ±0.013 

13.631 ± 0.046 

0.481 ±0.048 

-0.104 ±0.048 

0.584 ±0.018 

TrES-4B 

2013 

Jul 04 

K s 


13.5341 ± 0.0084 

13.777 ±0.024 


-0.243 ±0.025 


WASP-1B 

2012 

Jul 27 

K' 



14.827 ±0.077 




WASP-1B 

2013 

Aug 19 

K s 


15.130 ±0.046 

15.116 ± 0.055 


0.014 ±0.072 


WASP-1B 

2014 

Jul 12 

K s 



15.134 ± 0.064 




WASP-1B 

2014 

Oct 03 

K s 



14.837 ±0.100 




WASP-2B 

2012 

Jul 27 

K' 

13.046 ± 0.010 

12.4863 ±0.0061 

12.503 ± 0.028 

0.543 ±0.030 

-0.016 ±0.028 

0.560 ±0.012 

WASP-2B 

2013 

Jun 22 

K s 


12.4774 ± 0.0075 

12.319 ± 0.017 


0.158 ±0.018 


WASP-3B 

2012 

Jun 05 

K' 



15.889 ± 0.049 




WASP-3B 

2012 

Jul 27 

K' 

16.88 ±0.12 

16.090 ± 0.079 

16.002 ± 0.052 

0.87 ±0.13 

0.088 ± 0.095 

0.78 ±0.15 

WASP-3B 

2013 

May 31 

K s 


15.957 ±0.091 

15.913 ± 0.027 


0.044 ± 0.095 


WASP-8B 

2012 

Jul 27 

K' 

12.993 ± 0.062 

11.354 ±0.056 

10.636 ± 0.021 

2.357 ±0.065 

0.719 ±0.060 

1.638 ±0.083 

WASP-8B 

2013 

Aug 19 

K s 



10.646 ±0.018 




WASP-12B 

2012 

Feb 02 

K' 

14.188 ± 0.076 


13.51 ±0.11 

0.68 ±0.13 



WASP-12C 

2012 

Feb 02 

K' 

14.149 ±0.021 


13.765 ± 0.036 

0.384 ± 0.042 



WASP-12B 

2013 

Mar 02 

K s 



13.474 ± 0.030 




WASP-12C 

2013 

Mar 02 

K s 



13.370 ± 0.017 




WASP-14B 

2012 

Jul 27 

K' 

14.076 ± 0.020 

13.438 ± 0.024 

13.605 ± 0.025 

0.471 ± 0.032 

-0.167 ±0.034 

0.638 ±0.031 

WASP-14B 

2013 

Mar 26 

K' 



13.386 ± 0.057 




WASP-14B 

2013 

Mar 26 

K s 



13.322 ± 0.054 





References. — Primary star apparent magnitudes from 2MASS ijSkrutskie et aT1l2006T) . 

Note. — The K column indicates whether the following columns refer to the K' or K s bandpass, as both were used in the Friends of Hot 
Jupiters survey. The rrix columns refer to the secondary star’s apparent magnitude in the X filter, which is computed from the primary star’s 
apparent magnitude and the measured flux ratio as reported in Tabic [7T| The last three columns show the computed color of the companion star. 
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Table 5 

K-band astrometric measurements of confirmed common proper motion 

companions 


Companion 

UT Obs. Date 

K 

P 

(mas) 

PA 

(°) 

HAT-P-7B 

2012 Jul 27 

K' 

3858.7 ± 1.8 

89.958 ± 0.025 

HAT-P-7B 

2013 May 31 

K s 

3859.3 ± 1.7 

90.020 ±0.025 

HAT-P-7B 

2013 Jun 22 

K s 

3858.7 ±2.0 

90.020 ± 0.024 

HAT-P-7B 

2014 Jul 12 

K s 

3858.6 ± 1.7 

90.010 ±0.024 

HAT-P-8B 

2012 Jul 27 

K' 

1037.9 ± 1.5 

137.601 ± 0.084 

HAT-P-8C 

2012 Jul 27 

K' 

1047.8 ± 1.6 

140.868 ± 0.089 

HAT-P-8B 

2013 Aug 19 

K' 

1040.3 ± 1.6 

137.810 ± 0.084 

HAT-P-8C 

2013 Aug 19 

K' 

1045.0 ± 1.5 

141.058 ±0.083 

HAT-P-8B 

2013 Aug 19 

K s 

1041.0 ± 1.5 

137.770 ± 0.084 

HAT-P-8C 

2013 Aug 19 

K s 

1044.7 ± 1.5 

141.059 ±0.092 

HAT-P-10B 

2012 Feb 02 

K' 

342.5 ± 1.5 

214.09 ±0.24 

HAT-P-10B 

2013 Aug 19 

K' 

355.0 ± 1.5 

215.65 ±0.24 

HAT-P-10B 

2013 Aug 19 

K s 

355.0 ± 1.5 

215.64 ±0.24 

HAT-P-14B 

2012 Jun 05 

K' 

857.6 ± 1.5 

264.10 ±0.10 

HAT-P-14B 

2013 Mar 26 

K s 

857.4 ± 1.5 

264.24 ±0.11 

HAT-P-14B 

2014 Jul 07 

K s 

856.9 ± 1.6 

264.38 ±0.13 

HAT-P-16B 

2012 Feb 02 

K' 

691.6 ± 1.5 

153.83 ±0.12 

HAT-P-16B 

2013 Aug 19 

K s 

688.6 ± 1.5 

153.65 ±0.12 

HAT-P-24B 

2012 Feb 02 

K' 

4946.1 ± 1.8 

170.894 ± 0.020 

HAT-P-24B 

2014 May 13 

K s 

4944.2 ± 1.9 

170.872 ± 0.020 

HAT-P-30B 

2012 Feb 02 

K' 

3835.1 ± 1.7 

4.219 ±0.024 

HAT-P-30B 

2014 May 13 

Ks 

3836.6 ± 1.7 

4.206 ±0.024 

HAT-P-32B 

2012 Feb 02 

K' 

2935.5 ± 1.6 

110.583 ±0.030 

HAT-P-32B 

2013 Mar 02 

K s 

2935.5 ± 1.7 

110.624 ±0.031 

HAT-P-33B 

2012 Feb 02 

K' 

307.2 ± 1.5 

117.86 ±0.29 

HAT-P-33B 

2013 Mar 02 

K s 

306.3 ± 1.5 

118.05 ±0.29 

TrES-2B 

2012 Jul 27 

K' 

1105.4 ± 1.5 

136.325 ± 0.077 

TrES-2B 

2013 May 31 

Ks 

1106.7 ± 1.5 

136.357 ± 0.077 

TrES-4B 

2012 Jul 27 

K' 

1562.3 ± 1.6 

0.357 ± 0.056 

TrES-4B 

2013 Jul 04 

K s 

1563.4 ± 1.5 

0.386 ±0.055 

WASP-1B 

2012 Jul 27 

K' 

4580.0 ± 1.8 

1.901 ±0.021 

WASP-1B 

2013 Aug 19 

K s 

4581.5 ± 1.8 

1.953 ±0.021 

WASP-1B 

2014 Jul 12 

K s 

4582.3 ± 1.8 

1.958 ±0.021 

WASP-1B 

2014 Oct 03 

K s 

4581.7 ± 1.8 

1.917 ±0.021 

WASP-2B 

2012 Jul 27 

K' 

734.3 ± 1.5 

104.39 ±0.12 

WASP-2B 

2013 Jun 22 

K s 

730.0 ± 1.5 

104.54 ±0.12 

WASP-3B 

2012 Jun 05 

K' 

1192.2 ± 1.6 

87.103 ±0.077 

WASP-3B 

2012 Jul 27 

K' 

1191.0 ± 1.5 

87.070 ± 0.075 

WASP-3B 

2013 May 31 

K s 

1189.1 ± 1.6 

87.169 ±0.081 

WASP-8B 

2012 Jul 27 

K' 

4505.2 ± 1.7 

170.948 ± 0.021 

WASP-8B 

2013 Aug 19 

K s 

4507.4 ± 2.0 

170.988 ±0.021 

WASP-12B 

2012 Feb 02 

K' 

1058.8 ± 1.5 

251.242 ±0.084 

WASP-12C 

2012 Feb 02 

K' 

1067.1 ± 1.5 

246.751 ± 0.083 

WASP-12B 

2013 Mar 02 

Ks 

1058.6 ± 1.5 

251.444 ±0.082 

WASP-12C 

2013 Mar 02 

Ks 

1068.4 ± 1.5 

246.945 ± 0.082 

WASP-14B 

2012 Jul 27 

K' 

1450.8 ± 1.5 

102.066 ±0.059 

WASP-14B 

2013 Mar 26 

K' 

1449.1 ± 1.5 

102.210 ±0.060 

WASP-14B 

2013 Mar 26 

K s 

1449.5 ± 1.5 

102.207 ± 0.060 


Note. — We use both the K' and K s band for astrometric measurements. We imaged a few of our targets on the same night in both bands 
(HAT-P-8, HAT-P-10, and WASP-14) and show that the astrometry agrees to well within our uncertainties. The p column shows the on-sky 
separation between the fitted positions of the primary and companion stars, in milliarcseconds. The PA column shows the position angle, in 
degrees east of north. 
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Table 6 

Astrometric measurements from previous studies for confirmed common 
proper motion companions 





PA 


Companion 

U'l Obs. Date 

(mas) 

(°) 

Reference 


HAT-P-7B 

2009 Jul 18-22 a 

3870 ± 30 

90.4 ±0.5 

HAT-P-7B 

2009 Aug 06 

3875 ± 5 

89.8 ±0.3 

HAT-P-7B 

2009 Oct 30 

3820 ± 10 

90.4 ±0.1 

HAT-P-7B 

2009 Nov 02 

3861 ± 6 

89.8 ±0.1 

HAT-P-7B 

2011 Aug 12 

3871 ± 6 

89.7 ±0.1 

HAT-P-7B 

2011 Nov 09 

3860 ± 70 

89.9 ±0.8 

HAT-P-7B 

2012 Jul 07 

3860 ± 4 

89.8 ±0.1 

HAT-P-8B 

2012 Jun 24 

1040 ± 14 

137.9 ±0.8 

HAT-P-8B 

2013 Jul 03 

1053 ± 14 

137.6 ±0.8 

HAT-P-8C 

2012 Jun 24 

1049 ± 14 

141.4 ±0.8 

HAT-P-8C 

2013 Jul 03 

1041 ± 14 

140.7 ±0.8 

TrES-2B 

2007 May 3, 

1089 ± 8 

135.5 ±0.1 

TrES-2B 

2009 Jul 18-22 a 

1110 ±30 

137.0 ±2.0 

TrES-2B 

2009 Oct 29 

1085 ± 6 

136.1 ±0.2 

TrES-4B 

2008 Jun a 

1555 ± 5 

-0.2 ±0.1 

TrES-4B 

2009 Jul 18-22 a 

1540 ± 30 

1.2 ± 1.2 

TrES-4B 

2009 Oct 30 

1550 ± 7 

-0.1 ±0.2 

WASP-2B 

2007 Nov a 

757 ± 1 

104.7 ±0.3 

WASP-2B 

2009 Apr 13 

761 ±9 

103.5 ±0.2 

WASP-2B 

2009 Oct 29 

739 ± 24 

104 ± 1.3 

WASP-2B 

2011 Nov 09 

744 ± 13 

104.6 ±0.7 

WASP-12BC b 

2009 Oct 30 

1047 ± 21 

249.7 ±0.8 

WASP-12BC b 

2011 Nov 8-9 a 

1043 ± 14 

249.9 ±0.5 

WASP-12BC b 

2011 Dec 04 

1055 ± 26 

250.0 ± 1.0 

WASP-12BC b 

2012 Feb 25 

1078 ± 33 

249.4 ± 1.1 


^aed^et^L (2013) 
Narita et al. (2010b ) 
Bergfors et al (2013 ) 
Narit^e^aL ~ 
Narita^t^aL _ 
JBergfors^t^aL 

NaritaetaL (' 

BechteretaL 

BechteretaL 

BechteretaL 

BechteretaL 

Daemgen et al. 
^Faed^et^L (2013) 
Bergfors et al. (2013) 
Daemgen et al. (2009) 
^aedtet^U (2013) 

^er^for^t^L (JMS) 

^aemger^et^L (2009) 
^Bergfors^tiiii (2013) 
^ergfor^t^L (2013) 
< Ber^fors_et B ^L (2013) 
i Bergfors i ^t BB ^L (2013) 
^ergfor^t^L (2013) 

^rossfiel^^t_^L (2012) 

Crossfield et al. (2012) 


Note. — Astrometric measurements of Friends of Hot Jupiters companions reported in previous studies. Unless otherwise stated, only studies 
reporting a numerical value with uncertainties for both separation and position angle of all companions are included in this table and plotted in 
Figure [2] 

a A specific date was not reported. We report the date to the same precision as it appears in original work. 

k This study did not resolve the two companions. Thus, the measurement for the center of light of both companions is reported here and only 
used in the plots for the combined light of WASP-12BC. 


Table 7 

Astrometric measurements from previous studies for candidate companions 
determined to be background objects 


UT Obs. Date 

P 

(mas) 

PA 

(°) 

Reference 

For the second candidate companion to HAT-P-7 


2009 Aug 06 

3139 ± 5 

266.30 ± 0.37 

Narita et al. ( 

'2010b) 

2009 Nov 02 

3137 ± 6 

266.14 ±0.09 

Narita et al. ( 

2012) 

2011 Aug 12 

3103 ± 11 

266.23 ±0.24 

Narita et al. ( 

’20T2) 

2012 Jul 07 

3095 ± 4 

265.73 ± 0.08 

Narita et al. ( 

, 20T7) 

2012 Jul 27 

3091.5 ±3.3 

266.00 ±0.15 

this work 


2013 May 31 

3083.0 ± 2.7 

265.947 ± 0.059 

this work 


2014 Jul 12 

3065.5 ± 1.7 

265.705 ±0.044 

this work 


For the candidate companion to HAT-P-15 



2012 Feb 02 

7099.3 ± 9.9 

58.834 ± 0.097 

this work 


2014 Oct 03 

6382.6 ± 7.70 

58.193 ± 0.058 

this work 



Note, 


These measurements are used to generate the plots in Figure [3] 
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Table 8 

Derived stellar parameters of confirmed common proper motion companions 


Companion 

UT Obs. Date 


M 

log g 

D 

J-band T e ^ 

H-band T e ^ 

K-band T e ^ 

(K) 

(Mg) 

(cgs) 

(AU) 

(K) 

(K) 

(K) 

HAT-P-7B 

2012 Jul 27 

3321 ± 28 

0.215 ±0.016 

5.047 ±0.018 

1238 ± 67 

329ll“ 

3357li| 

3316.71“ 

HAT-P-7B 

2013 May 31 

3307 ±10 

0.2089 ± 0.0043 

5.0563 ± 0.0047 

1238 ± 67 

3307.0 ±6.5 

HAT-P-7B 

2013 Jun 22 

3320 ± 30 

0.214 ±0.018 

5.047 ±0.019 

1238 ± 67 


335ll“ 

3291 I 2 ® 

HAT-P-7B 

2014 Jul 12 

3328 ± 11 

0.2215 ±0.0066 

5.0428 ± 0.0070 

1238 ± 67 


3328.11“ 

HAT-P-8B 

2012 Jul 27 

3216 ± 59 

0.168 ±0.020 

5.102 ± 0.029 

236 ± 13 

3153 

3293.2l° 1 8 

32081“ 

HAT-P-8B 

2013 Aug 19 

3252 ± 41 

0.182 ±0.017 

5.086 ±0.022 

237 ± 13 


3306.2lg'g 

3216l 2 g 

HAT-P-8C 

2012 Jul 27 

3058 ± 63 

0.131 ±0.011 

5.167 ± 0.023 

239 ± 13 

2989i4° 

3136ljg 

30551“ 

HAT-P-8C 

2013 Aug 19 

3130 ± 44 

0.146 ±0.010 

5.139 ±0.018 

238 ± 13 


31881 jg 

30921“ 

HAT-P-10B 

2012 Feb 02 

3483 ± 43 

0.344 ± 0.047 

4.922 ± 0.038 

41.7 ± 1.4 

3525.41®;® 

3443.llj g 

HAT-P-10B 

2013 Aug 19 

3494 ± 37 

0.362 ±0.035 

4.913 ±0.031 

43.2 ± 1.5 

3545.Ol^J 

3469. 2 I ® 2 

HAT-P-14B 

2012 Jun 05 

3310 ± 17 

0.211 ±0.010 

5.054 ±0.011 

175.8 ±9.4 



3310.2l®“ 

HAT-P-14B 

2013 Mar 26 

3356 ± 57 

0.232 ±0.037 

5.021 ± 0.039 

175.8 ±9.4 


3413lig 

33021“ 

HAT-P-14B 

2014 Jul 07 

3263 ± 20 

0.1873 ±0.0087 

5.082 ± 0.011 

175.7 ±9.4 


3263.llgg 

HAT-P-16B 

2012 Feb 02 

3255 ±67 

0.182 ±0.027 

5.081 ± 0.035 

162.5 ± 6.9 

3319.2 ± 7.3 


3194.81®“ 

HAT-P-16B 

2013 Aug 19 

3269 ± 51 

0.188 ±0.023 

5.076 ± 0.028 

161.8 ±6.9 


3315.11^ 

3225 ± 15 

HAT-P-24B 

2012 Feb 02 

3434 ± 26 

0.298 ±0.026 

4.968 ± 0.022 

1959 ± 99 



34341“ 

HAT-P-24B 

2014 May 13 

3499 ± 29 

0.373 ± 0.030 

4.908 ± 0.025 

1958 ± 99 

3474.41“ 

3498.8l|^ 

3523.2 ± 1.6 

HAT-P-30B 

2012 Feb 02 

3634 ± 29 

0.489 ± 0.020 

4.808 ± 0.020 

740 ± 31 


36341“ 

HAT-P-30B 

2014 May 13 

3692 ± 40 

0.519 ±0.019 

4.777 ± 0.019 

740 ± 31 

3709.3 ± 2,4 

3692.2 ± 1.6 

36741“ 

HAT-P-32B 

2012 Feb 02 

3516 ± 12 

0.393 ±0.012 

4.8930 ± 0.0098 

831 ± 15 

3524 . 0 I 49 


3508.6 ±6.5 

HAT-P-32B 

2013 Mar 02 

3551 ± 10 

0.4243 ± 0.0085 

4.8677 ± 0.0070 

831 ± 15 

3555li4 

35481“ 

HAT-P-33B 

2012 Feb 02 

3653 ± 54 

0.493 ± 0.031 

4.800 ± 0.031 

118.9 ±2.8 

3705ljg 


3604l 2 ° 

HAT-P-33B 

2013 Mar 02 

3776 ± 32 

0.557 ± 0.014 

4.740 ±0.013 

118.5 ±2.8 

38001“ 

37531“ 

TrES-2B 

2012 Jul 27 

3669 ± 29 

0.509 ±0.013 

4.787 ± 0.014 

216 ± 13 

3651.6llj 

3701.2 ± 2.4 

3655.71°® 

TrES-2B 

2013 May 31 

3662 ±39 

0.501 ±0.021 

4.793 ± 0.021 

216 ± 13 

3697.91^1 

3627.2 ±8.1 

TrES-4B 

2012 Jul 27 

3985 ±93 

0.626 ±0.028 

4.677 ± 0.023 

900 ± 89 

4048.2 ± 6.5 

4023.8 ±6.5 

38891“ 

TrES-4B 

2013 Jul 04 

3910 ± 110 

0.602 ± 0.036 

4.696 ± 0.032 

900 ± 89 


4 OI 9 . 8 II 3 

3817±11 

WASP-1B 

2012 Jul 27 

3446 ± 15 

0.313 ±0.019 

4.956 ±0.014 

1587±^|° 


3446.01“ 

WASP-1B 

2013 Aug 19 

3420 ± 26 

0.285 ±0.023 

4.978 ±0.021 

15871“° 


3441.91®;® 

3398.0 ±8.1 

WASP-1B 

2014 Jul 12 

3395 ±17 

0.266 ±0.014 

4.998 ±0.013 

1587±^|° 


3394.8l“ g 

WASP-1B 

2014 Oct 03 

3442 ± 18 

0.308 ±0.021 

4.960 ±0.017 

15871:“° 



34421“ 

WASP-2B 

2012 Jul 27 

3509 ± 26 

0.383 ± 0.028 

4.899 ± 0.022 

113.0 ±6.1 

3527.9 ± 2.2 

3523.55114 2 

3477.4l^J 

WASP-2B 

2013 Jun 22 

3514 ± 16 

0.391 ±0.016 

4.895 ± 0.013 

112.3 ±6.1 


3525.7l^J 

3502.6 ±2.9 

WASP-3B 

2012 Jun 05 

2909 ± 33 

0.1089 ±0.0035 

5.223 ±0.010 

300 ± 22 



29091^ 

WASP-3B 

2012 Jul 27 

2871 ± 51 

0.1053 ±0.0048 

5.233 ±0.015 

299 ± 22 

2825 ± 47 

2924l 2 ^ 

2867l 2 g 

WASP-3B 

2013 May 31 

2922 ± 48 

0.1109 ±0.0058 

5.216 ±0.016 

299 ± 22 


29621 2 ® 

2883.6 ±9.8 

WASP-8B 

2012 Jul 27 

3530 ± 130 

0.34 ±0.11 

4.881 ± 0.099 

383 ± 48 

3383.3 ±9.6 

3568 ± 12 

3666.0l®2 

WASP-8B 

2013 Aug 19 

3657 ± 76 

0.504 ± 0.043 

4.793 ±0.043 

384 ± 48 



3657.2l|| 

WASP-12B 

2012 Feb 02 

3786 ± 53 

0.561 ±0.022 

4.736 ± 0.021 

462 ± 39 

3810l^| 


3762 I 35 

WASP-12B 

2013 Mar 02 

3769 ± 44 

0.554 ±0.020 

4.743 ±0.018 

462 ± 39 



3768.6l“g 

WASP-12C 

2012 Feb 02 

3748 ± 90 

0.540 ± 0.038 

4.753 ± 0.037 

466 ± 40 

3828.0 ±9.8 


36741“ 

WASP-12C 

2013 Mar 02 

3808 ± 49 

0.571 ±0.019 

4.727 ± 0.017 

466 ± 40 



3807.6lgg 

WASP-14B 

2012 Jul 27 

3464 ± 31 

0.331 ± 0.037 

4.939 ± 0.028 

302 ± 23 

3476.0 ± 3.3 

3484.211^ 

3432.213'^ 

WASP-14B 

2013 Mar 26 

3464 ± 19 

0.336 ± 0.024 

4.939 ±0.018 

302 ± 23 



3461.4l 8 ;g 


Note. — The T e ff, M, log p, and D columns are error weighted averages from all measurements of each target for each date. These uncertainties 
include measurement uncertainty, uncertainties from the primary star’s stellar parameter and uncertainties on the error weighted average but do 
not include uncertainties arising from the stellar models and our assumptions on stellar composition. The uncertainties in last three columns, 
which report T e ff for each filter, only include measurement uncertainties and thus would be underestimates of the true uncertainty. The last three 
columns show that the companion object have T e ff measurements consistent with a late type main sequence star in all filters. 
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Table 9 

Stellar companion fractions from other surveys 


Survey 

V 

Range 

^targets 

Notes 

This work 





Total sample, uncorrected 

34% ± 7% 

5.5" 

50 

Diffraction-limited AO 

Total sample, corrected 

49% ± 9% 

5.5" 

50 

Diffraction-limited AO 

Misaligned sample, corrected 

48% ± 12% 

5.5" 

27 

Diffraction-limited AO 

Control sample, corrected 

51% ± 13% 

5.5" 

23 

Diffraction-limited AO 


Near-infrared direct 


imaging surveys of transiting gas giant planet hosts 



21 % ± 12% 

Unknown 3, 

14 

Lucky imaging 

38% ± 15% 

Unknown 3, 

16 

Lucky imaging 

29% ± 12% 

Unknown 3, 

21 

Lucky imaging 

33% ± 15% 

Unknown 3, 

15 

Diffraction-limited AO 


Near-infrared direct imaging surveys of Kepler planet 


candidate hosts 


( Adams >i etjaL (2012) 59% d= 5% 

Adams et al. (20135 33% d= 17% 

Dressing et al7 72014) 31% d= 5% 

Lillo-Box et al. (2012 . 2014) 33% d= 4% 

Wang et al"72014b) 45% d= 7% 


6 " 

90 

Diffraction-limited AO 

Unknown 3, 

12 

Diffraction-limited AO 

4" 

87 

Diffraction-limited AO 

6 " 

174 

Lucky imaging 

6 " 

56 

Diffraction-limited AO 


Optical direct imaging surveys of Kepler planet candidate hosts 


Law et al. (2014) 

7.4% ± 1.0% 

2.5" 

715 

Diffraction-limited AO 


14,0 ±3 5% 



Only hosts of planets P < 15 days, R p > -RNeptune 

Gilliland et al. (2015) 

22 % ± 10% 

2.5" 

23 

Hubble Space Telescope WFC3 imaging 

Horch et al. (2014). raw 

7.0% ± 1.1% 

1 " 

588 

Differential speckle survey on WIYN 3.5m telescope 


23% ± 8% 

1 " 

35 

Differential speckle survey on Gemini-N 8.1m telescope 

Horch et al. (2014). corrected 

37% ± 7% 

1 " 

588 

Differential speckle survey on W1YJN 3.5m telescope 


47% ± 19% 

1 " 

35 

Differential speckle survey on Gemini-N 8.1m telescope 

Multi-modal surveys of field stars 




Raghavan et al. (2010) 

44% ± 2% 
24% ± 1% 

varies* 5 

See note c 

454 

454 

Solar-like stars with 25 pc 

For comparison with our survey 


Note. — Unless otherwise noted, companion fraction rj reported are based on raw counts, with no correction for completeness. The uncertainties 
for completeness corrected companion fractions do not scale inversely with the square root of the number of targets because the completeness 
correction effectively multiplies both the value and uncertainty by a constant. 

a These studies do not report a number of companions in any given range. 

k This study uses multiple techniques to find stellar companions, each with different sensitivities for companion separation. 

C This line reports the fraction of stars with companions in the same range as our survey, as detailed in Section l5.3l 




















